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ABSTRACT 


PEAPS  (Passive  Environmental  ASW  Prediction  System)  is 
a relatively  unsophisticated  model  which  accepts  input  source 
and  receiver  parameters  and  then  predicts  sound  propagation 
characteristics  in  an  ocean  environment,  the  corresponding 
transmission  loss,  and  the  probability  of  detection.  The 
program  was  written  for  a programmable  desk-top  calculator 
for  immediate  deployment  and  operational  testing  aboard 
small  ASW  platforms.  The  program  is  also  available  in  a 
form  suitable  for  larger  computers. 
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CD 3 to  CD 5 
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CZSTOP 

CZSTRT 

DFL 


DEFINITION 

The  absorption  coefficient  (dB/yd) 
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The  angle  of  incidence  below  which  the 
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G 
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DEFINITION 
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for  follow-on  programs. 

R 

The  horizontal  range  of  the  ray  that 
is  being  traced  (in  ft)  and  the  TL/PD 
range  (in  yards) . 

RD 

The  input  receiver  depth  (in  feet) . 

RADINV 

The  inverse  of  the  radius  of  curvature 
for  a ray  as  it  is  traced  through  a 
vertical  section  of  the  water  column. 
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SYMBOL 


DEFINITION 


KPT 

RS 

RT 


SD 

SE 

SF 

SS 

SSL 

STL 


SV 


The  transition  range  used  in  the  deep- 
sound-channel  segment  of  the  program. 

The  derived  skip  distance  used  in  the 
surface-duct  segment  of  the  program. 

The  derived  transition  range  used  in 
the  surface-duct  segment  of  the 
program. 

The  input  source  depth,  in  feet. 

The  derived  signal  excess 
(SE  = STL  - IFOM) . 

The  input  source  frequency,  in 
Hertz  (Hz) . 

Speed  of  sound. 

The  derived  sea  surface  loss  used  in 
the  surface  duct  segment  of  the  program. 

The  cumulative  transmission  loss  from 
the  existing  transmission  paths  at  a 
particular  range  from  the  source. 

The  derived  speed  of  sound  at  the 
source. 


SI  to  S15  The  derived  speed  of  sound  corresponding 

to  the  input  data  from  the  BT  trace. 

TH  The  angle  the  ray  makes  with  the 

horizontal  as  it  is  being  traced. 

THBB  The  derived  minimum  emission  angle 

from  the  source  that  will  result  in 
a bottom  bounce  ray. 

THBL  The  angle  used  to  determine  the  bottom 

loss  in  the  bottom  bounce  segment  of 
the  program. 

THMAX  The  maximum  ray  emission  angle  from 

the  source  entered  by  the  operator 
when  MODE  = 2 . 

THMIN  The  minimum  ray  emission  angle  from 

the  source  entered  by  the  operator 
when  MODE  = 2. 


svMRrvr, 


DEFINITION 


THP 


THPP 

TL 

WRKNG 

Z 


The  approximate  minimum  source 
emission  angle  that  will  produce  a 
convergence  zone  ray.  Used  as  a 
starting  point  in  the  automatic  modes 
(MODE  = 0 or  1)  to  determine  CZMIN 
by  iterative  ray  tracing. 

The  approximate  THBB  used  in  the  auto- 
matic modes  as  a starting  point  for 
iterative  ray  tracing  to  determine 
THBB. 

A transmission  loss  storage  location. 

Storage  registers  used  in  several 
portions  of  the  programs  for 
temporary  storage. 

The  depth  of  the  ray  being  traced, 
measured  from  the  source  depth. 
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I . INTRODUCTION 


A.  BACKGROUND 

Shore-based  prediction  systems  rely  on  complex  acoustic 
models  which  require  large  computers  and  relatively  long 
computation  times.  Their  predictions  are  based  on  histori- 
cal meteorologic  and  oceanographic  data  for  regions  in  which 
environmental  conditions  are  assumed  to  be  homogeneous . 

These  systems  use  a single  level  of  sonar  performance  as  the 
basis  for  all  range  predictions.  The  value  of  such  systems 
lies  in  their  use  as  research  tools  and  for  pre-deployment 
planning:  they  can  be  considered  to  be  strategic  systems. 

With  the  advent  of  the  programmable  desk-top  calculator 
comes  the  possibility  that  a simple  acoustic  prediction  model 
can  be  placed  on  board  anti-submarine  warfare  (ASW)  platforms. 
This  makes  it  possible  for  the  ASW  platofrm  to  take  advantage 
of  on-scene  bathythermographic  (BT)  data  to  obtain  current 
range  predictions  based  on  the  output  of  the  platform's  own 
sensors.  In  addition,  an  on-board  system  enables  the  ASW 
platform  to  obtain  predictions  under  all  electronic  emission 
control  (EMCON)  conditions.  The  programmable  desk-top  cal- 
culators make  possible  environmental  predictions  at  the 
tactical  level. 
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B.  PEAPS 


1.  Description 

PEAPS  (Passive  Environmental  ASW  Prediction  System) 
was  designed  for  use  in  any  programmable  desk-top  calculator 
or  general-purpose  computer  with  sufficient  capacity  for 
the  program. 

The  system  is  composed  of  four  sub-programs:  a 

speed  of  sound  profile  (SSP)  generator,  a ray  trace  generator, 
a transmission  loss  (TL)  generator,  and  a probability  of 
detection  (PD)  generator.  The  SSP  and  ray  trace  generators 
are  combined  as  an  SSP/Ray  Trace  program;  the  TL  and  PD 
generators  are  combined  as  a TL/PD  program.  These  sub- 
programs are  described  in  detail  in  subsequent  chapters. 

2 . Assumptions  and  Constraints 

The  system  was  written  originally  in  a machine 
language  for  use  on  a small  capacity,  programmable,  desk- 
top calculator  (the  Hewlett-Packard  Model  9810A  with  expanded 
storage) . The  hardware  used  to  develop  PEAPS  is  described  in 
Appendix  A.  Total  program  capacity  of  the  calculator  was 
2036  program  steps  with  100  data-storage  registers.  This 
calculator  was  used  because  it  was  relatively  inexpensive 
and  readily  accessible;  it  was  reasoned  that  if  a simple 
working  model  could  be  designed  on  a small-capacity, 
programmable  calculator,  more  sophisticated  models  would  be 
possible  on  larger  programmable  calculators. 

PEAPS  is  a relative  unsophisticated  system  which 
uses  approximation  techniques  to  compute  transmission  losses 
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under  the  assumptions  of  a flat-bottomed  ocean  and  phase 
incoherent  interference.  Lateral  variations  in  vertical 
temperature  profiles  are  neglected.  Speed  of  sound  profiles 
are  approximated  by  straight-line  segments,  with  storage 
available  for  fourteen  segments  (fifteen  input  points)  and 
five  smoothed  gradients.  When  available  to  the  user  deep 
ocean  historical  data  may  be  incorporated  into  the  input, 
but  such  data  are  neither  stored  nor  generated  internally. 

The  specific  assumptions  and  constraints  used  in  each  sub- 
program are  described  in  the  subsequent  chapter  which  deals 
with  that  sub-program. 

General  worksheets  for  PEAPS  have  been  designed  to 
guide  the  operator  through  the  sequence  of  steps  required 
to  generate  output  (Appendix  B) . 

Operator  note  sheets  also  have  been  written  in 
outline  form  to  describe  the  use  of  the  individual  sub-programs 
of  the  system  and  to  assist  the  operator  when  problems  and 
questions  arise  (Appendix  C) . 
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II.  SPEED  OF  SOUND  PROFILE  (SSP)  GENERATOR 


A.  PURPOSE 

The  speed-of-sound  generator  produces  a plot  of  sound 
speed  versus  depth  for  inputs  from  a bathythermograph  (BT) 
trace.  The  data  are  then  searched  for  critical  points  which 
are  stored  for  subsequent  use  by  the  ray  trace  and  transmis- 
sion loss  (TL)  programs.  A critical  point  is  a point  on  the 
sound  speed  profile  where  the  gradient  of  the  sound  speed 
changes  sign  or  becomes  zero.  All  sound-speed  points  are 
stored  for  use  in  the  ray  trace  sub-program. 

B.  INPUTS 

Figure  1 shows  the  SSP  portion  of  the  operator's  work- 
sheet. The  operator  manually  inserts  the  BT  data  from  the 
worksheet  into  the  calculator  programmed  for  SSP/Ray  Trace. 

As  indicated,  the  data  entered  on  the  worksheet  are  the 
number  of  BT  points  to  be  entered,  and  the  depth  (ft) , 
temperature  (°F) , and  salinity  (ppt)  of  each  input  point; 
if  salinity  is  unknown  and  historical  salinity  data  are 
unobtainable,  a value  of  35  parts  per  thousand  (ppt)  is 
entered.  If  the  bottom  depth  is  greater  than  that  of  the 
BT  trace  and  historical  data  are  not  available,  it  is  recom- 
mended that  the  operator  enter  the  last  BT  temperature  as 
the  bottom  temperature.  This  will  produce  a profile  for 
constant  temperature  from  the  last  BT  point  to  the  bottom. 

If  tables  of  historical  salinities  and  temper at ures-versus-depth 
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PEAPS  OPERATOR  WORKSHEET 


(See  PEAPS  Operator  Note  Sheets) 

BT  Information  DTG: Z 

Location:  N/S , E/W 

1.  Turn  the  calculator  and  plotter  on.  "LOAD"  the  SSP/Ray 
trace  program  cards.  Put  the  plotting  paper  on  the 
plotter  and  set  the  "GRAPH  LIMITS"  of  the  plotter. 

2.  Press  "STOP,  END,  CONTINUE"  (this  sequence  is  used  when 
a new  SSP  is  to  be  generated. 


Note  1:  If  you  enter  an  erroneous  value: 

a.  If  you  have  not  pressed  "CONTINUE",  press 
"CLR  X",  then  enter  the  correct  value. 

b.  If  you  have  pressed  "CONTINUE",  you  must 
re-enter  all  the  points  again;  press  "GO  TO, 
LABEL,  A" , then  continue  with  Step  3 below. 

3.  Press  "CONTINUE". 

4.  SSP  inputs:  ENTER  number  of  BT  points  (15  maximum) 

Press  "CONTINUE"  after  each  value  is  entered  to  advance 
the  program  output 

Input  Input  Input  Speed  of  Sound 

Point  Depth  (ft)  Temp  (°F)  Salinity  (ppt)  (ft/sec) 
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FIGURE  1 


for  the  oceanic  area  of  interest  are  available  to  the 
operator,  then  use  Section  G of  the  Operator  Note  Sheets 
(Appendix  G)  to  merge  the  historical  and.  on-scene  BT  data 
for  improved  accuracy. 

C.  MODEL 

The  basic  physical  concepts  of  sound  transmission  are 
fairly  simple,  but  computation  of  the  sound  speed  for  the 
ocean  is  rather  complex  because  of  the  varying  physical 
properties  of  sea  water.  The  computation  of  sound  speed  from 
empirical  relations  based  upon  laboratory  measurements  began 
in  the  early  1950's  (Del  Grasso,  1952;  Wilson,  1960) . 

Attempts  to  record  accurate  hydrographic  data  while  elec- 
tronically measuring  the  sound  speed  were  successful  in  the 
early  1960's  (Hays,  1961;  Mackenzie,  1961).  For  these  latter 
experiments,  a computed  sound  speed  was  obtained  based  on 
three  measured  parameters:  temperature,  salinity,  and  depth. 

The  results  were  then  compared  with  the  electronically- 
recorded  sound  speeds.  These  observations  resulted  in  either 
slight  modification  of  the  coefficients  in  existing  equations 
or  the  development  of  new  polynomial  equations  (Wilson,  1960; 
Leroy,  1968;  Frye  and  Pugh,  1971;  Anderson,  1971;  Del  Grasso, 
1972).  From  the  results  obtained  by  these  experimenters,  the 
variation  of  sound  speed  versus  temperature,  salinity,  and 
depth  can  be  approximately  represented  as  indicated  in  Table  I. 

Approximately  99%  of  all  sea  water  falls  within  the 
temperature  range  -3°C  to  +30  °C  and  the  salinity  range  33  to 
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TABLE  I 


Var!  it ions  of  sound  speed  caused  by  variations  in 
temperature,  salinity,  and  depth.  A variation  of  one 
parameter  signifies  that  the  other  parameters  are  held 
constant:  temperature  at  0°C,  zero  depth,  salinity  at 

35.00  ppt. 

Temperature  (°C) 

0-10 
10  - 20 
20  - 30 

Salinity  (ppt) 

33-35  1.55  (meters/sec) /ppt 

35  - 37  1.25  (meters/sec/ /ppt 


Sound  Speed  Change 
4.5  (meters/sec) /°C 

3.0  (meters/sec) /°C 

2.0  (meters/sec) /°C 


Depth  (meters) 

0 - 10,000 


0.017  (meters/sec) /meter 
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37  ppt  (Montgomery,  1958;  Frye  and  Pugh,  1971).  Thus, 
thermal  gradients  are  usually  much  more  significant  than 
salinity  gradients  in  sound  speed  computations;  so  that  the 
temperature-depth  relation  is  a factor  which  must  be 
considered  in  greatest  detail  in  any  sound  propagation  model. 

Because  the  Frye  and  Pugh  equation  [1]  for  sound  speed 
is  dominated  by  temperature  terms  and  requires  relatively 
few  program  steps  compared  with  other  speed  of  sound  equations, 
it  was  used  in  the  sound  speed  profile  sub-program.  The 
equation  was  modified  to  use  the  approximation 


ASS  _ ASS  AP  .i  A i ASS 
AZ  AP  AZ  U,X  AP 


for  the  speed  of  sound  SS  in  meters/sec,  the  pressure  P in 
bars  (gage  pressure) , and  the  depth  Z in  meters.  The 
resulting  equation  is 


SS  = 1449.3  + 1.5848  x 10~2Z  + 1.572  x 10_7Z2 

- 3.46  x 1<T16Z4  + 4.587T  - 5.356  x 10"2T2 

+ 2.604  x 10"4T3  + 1 . 19  (S  - 35)  + 9.6  x 10~2 (S  - 35) 2 
+ 1.354  x 10~5T2Z  - 7.19  x 10~9TZ2 

- 1.2  x 10~2(S  - 35) T 


where  SS  = sound  speed  in  meters/sec 
Z = depth  in  meters 
T = temperature  in  °C 

S = salinity  in  ppt  (parts  per  thousand)  «. 
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The  SSP  sub-program  converts  the  input  data  to  MKS  units 
for  use  in  the  above  equation  and  then  converts  the 
resultant  sound  speeds  to  feet/sec  for  use  in  subsequent 
routines. 

The  above  equation  is  relatively  inaccurate  in  areas 
where  salinity  exceeds  36.6  ppt:  Frye  and  Pugh  claim 

accuracy  only  to  36.6  ppt  and  are  presently  working  to  extend 
the  domain  of  their  equation. 

D.  OUTPUT 

The  calculator  prints  the  input  data  and  the  corresponding 
sound  speeds  on  tape  (Fig.  2) , makes  a plot  of  sound  speed 
versus  depth  (Fig.  3) , and  refers  the  operator  to  Operator 
Note  2 for  ray  trace  preparation  (see  Appendix  C) . 
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figure  2.  Sample  Speed  of  Sound 
Profile  generator 
printer  input/output. 
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III.  RAY  TRACE  GENERATOR 


A.  PURPOSE 

A ray  trace  sub-program  included  in  PEAPS  makes  the 
following  tests  and  decisions: 

1.  inspect  the  speed-of-sound  profile  to  determine 
whether  the  characteristics  necessary  for  convergence  zone 
propagation  are  present  or  not; 

2.  if  CZ  conditions  exist,  determine  the  maximum  and 
minimum  CZ  ranges,  where  minimum  CZ  range  is  defined  as  the 
range  from  the  source  to  the  inner  limit  of  the  CZ  annulus 
and  maximum  CZ  range  is  defined  as  the  range  to  the  outer 
limit  of  the  CZ  annulus  (Pig.  4),  and  in  any  event; 

3.  determine  the  minimum  bottom  bounce  angle. 

B.  MODES 

The  ray  trace  sub-program  has  three  modes  of  operation 
designated  by  0,  1,  and  2. 

Mode  0 automatically  determines  and  plots  the  limiting 
rays  of  the  convergence- zone  ray-family  and  stores  the 
corresponding  ranges. 

Mode  1 is  the  same  as  Mode  0,  but  omits  the  ray  trace 
plot. 

Mode  2 gives  a plot  of  the  operator's  input  ray-family 
followed  by  a Mode  1 run. 
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Convergence  Zone  Path 


C.  INPUTS 


The  operator  enters  the  Mode  number  and  source  depth 
from  the  ray  trace  portion  of  the  worksheet.  If  Mode  2 
is  selected,  he  must  also  enter  twice  the  maximum  and 
minimum  emission  angles  and  twice  the  desired  increment  in 
emission  angle  between  successive  rays1  (Fig.  5) . 

D.  MODEL 

The  ray  trace  generator  uses  Snell’s  Law  in  a small 
angle  approximation. 

Snell's  Law  for  rays  can  be  written  as 

Ar  = R(sin02  - sin81) 

and 

AZ  = R(cos01  - cos02) 

where  the  geometry  is  suggested  in  Fig.  6.  Elimination  of 
02  gives  the  range-depth  equation  in  terms  of  the  initial 
angle  of  depression, 

AZ1  = Rcos01  - R[1  - (sinOl  + ~)2}h  (1) 


1See  Section  B of  Note  2 in  the  Operator  Note  Sheets 

(Appendix  C) . 
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Note  2:  Change  the  plotting  paper  on  the  plotter.  Enter 

the  ray  trace  program  MODE  indicator. 

0 = Automatic  compute  and  plot 

MODE:  1 = Automatic  compute,  no  plot 

2 = Manual  inputs  for  plotting,  Automatic 
compute . 

If  a ray  trace  input  error  is  made,  press  "GOTO, 
LABEL,  B,  CONTINUE",  then  input  the  MODEi,  indicator 
and  ray  trace  data  again.  (This  sequence  is  used 
to  generate  a new  trace  using  the  old  SSP  aata) 

Press  "CONTINUE" 


4,  Ray  trace  inputs: 

Press  "CONTINUE"  after  each  input  to  advance  the  program. 

Source  Depth:  (ft) 

If  Mode  2 is  used: 


Minimum  Ray  Angle  x 2: 


(degrees) 


Angle  Increment  x 2: (degrees) 

Maximum  Ray  Angle  x 2:  (degrees) 


Note  3:  If,  while  in  Mode  0 and  after  the  bottom  limiting 

ray  has  been  drawn,  a channelled  ray  is  drawn  to 
the  plot's  range  limit;  see  Operator  Note  sheet 
section  Note  2E. 

Change  the  plotting  paper  on  the  plotter.  Press 
"STOP,  END".  LOAD  the  TLOSS/PD  program  cards. 
Press  "STOP, END".  (The  "STOP,  END"  sequence  is 
used  when  a TLOSS/PD  input  error  is  made,  or  a new 
TLOSS/PD  -lot  is  to  be  generated  based  on  the  old 
SSP/Ray  trace  data.) 


5.  TLOSS/PD  inputs: 

Press  "CONTINUE"  after  each  input  to  advance  the  program. 
Press  "CONTINUE" 

Source  Frequency:  (Hz) 

Figure  of  Merit  (FOM)  : (dB) 

System  Sigma:  (use  10  if  unknown) 

Wave  Height:  0 1 2 3 4 5 6 7 8 9 10  11  12  (ft) 

Receiver  Depth:  (ft) 

Bottom  Type:  123456789 


FIGURE  5 
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Expanding  Azi  as  a MacLaurin  series  about  Ar  = 0 gives 

AZ2  = tanei(Ar)  + =■  (Ar)2 

2R(cos61)  J 

+ 1.-1 isnml  (Ar)4  + . . 

8RJ (cos01)  ' 

Substitution  of  the  trigonometric  series  representations 
of  tan  x and  cos  x into  the  first  two  terms  of  equation  (2) 
gives 

AZ3  = 61(1  + $~)  (Ar)  + (1-+  '2R^1~)  {Ar)2  * (3) 

To  conserve  program  space  and  reduce  the  computation 
time,  PEAPS  uses  a simplified  form  of  Eq.  (3) . Use  of  the 
small  anglve  approximations 


sin01 


2R2 (cos  61) 5 


(Ar) 


(2) 


and 


sin8  = 0 

0* 

cose  = l - , 


where  0 is  in  radians,  yields  the  simple  approximations 


AZ 


(61)  (Ar)  + 


(Ar)2 

R 


(4) 


and 


A0  = Ar/R  . 
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For  Eq.  (4)  to  accurately  represent  Eq.  (3) , restrictions 
must  be  imposed  on  Ar  and  Gl.  In  order  that  omitting  01 
introduces  no  more  than  10%  error,  01  must  be  less  than 

0.33  radians.  To  determine  the  limits  on  Ar,  examine  the 
case  of  "worst  possible"  conditions: 

1.  A maximum  gradient  of  0.5  sec  ^ and  a minimum  sound 
speed  of  4700  ft/sec  are  assumed  so  that  the  minimum  radius 
of  curvature  is  9400  ft. 

2.  The  maximum  allowable  angle  of  depression  of  0.35 
radians  (20°)  is  taken. 

To  ensure  that  no  layer  is  missed  or  poorly  represented  in 
the  ray  trace,  five  depth  increments  are  required  within  any 
layer  of  depth  D.  Then,  the  first-order  term  in  Eq.  (4) 
yields  the  equation 


Ar  < 


D 

501 


0.57D 


(5) 


as  an  approximate  bound  on  the  allowed  Ar  within  a layer 
of  depth  D. 

So  that  the  truncation  of  Eq.  (2)  to  second  order  in  Ar 
introduces  no  greater  error,  the  third  term  of  the  equation 
must  be  much  less  than  the  second, 


^'--r  Ur)3  « 


2R  (cosol) 


2R(cos01) 


(Ar) 


or 


Ar  « --f ■ < 2.6R 
smei  ■» 
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If  Ar  is  required  to  be  at  least  one  order  of  magnitude 
smaller  than  2.6R,  then 

Ar  < 0.26R  <_  2444  ft. 

However,  since  for  small  01  the  third  term  becomes 
arbitrarily  small,  it  is  also  necessary  to  require  that  the 
fourth  term  of  Eq.  (2)  be  much  less  than  the  second, 

1 (Ar)4  « A (Ar)2  , 

8R3(cos01)/  2R(cos01) 3 


or 


Ar  « (2R)  (cos  61)  2[1  + 4(sin81)2]"5  = 2.1R  . 

If  Ar  is  to  be  at  least  one  order  of  magnitude  less 
than  2.1R,  then 

Ar  < 0.21R  < 2012  ft. 

To  estimate  the  error  in  calculating  AZ  from  Eq.  (4) , 
subtract  Eq.  (4)  from  Eq.  (3): 

AZ3  - AZ  = — (Ar)  + (Ar)  2 


(6) 


The  fractional  error  is  given  by: 


AZ3  - AZ  _ 0.014  + (9.56  x l(Tb) (Ar) 
AZ3  0.36  + (1.16  x 10"4) (Ar) 


Away  from  a turning  point,  for  large  angles  where  the 
difference  defined  by  Eg.  (6)  is  large,  the  error  is 
approximately  4%  for  small  Ar  and  approximately  8%  for  large 
Ar.  If  the  error  is  required  to  be  no  greater  than  5%,  then 
Ar  must  be  no  greater  than  about  1000  feet  (near  a turning 
point,  AZ3  - AZ  approaches  zero  very  rapidly,  so  that  errors 
in  depth  are  negligible  for  small  angles  of  elevation  or 
depression) . 

Based  on  these  error  requirements,  the  simple  equation 


Ar  = 5(D) 


was  used  for  Ar  and  D in  feet.  For  the  ranges  of  D employed 
in  PEAPS  the  restrictions  on  Ar  set  forth  above  are  not 
seriously  violated  (Table  II) . 


E . OUTPUT 

If  convergence  zone  (CZ)  conditions  are  not  present, 
"NO  CZ"  appears  on  the  printer.  In  all  cases  "SEE  NOTE  3" 
and  "END"  appear,  indicating  the  end  of  the  SSP/Ray  Trace 
programming  sequence,  and  referring  the  user  to  Operator 
Note  3 (Appendix  C) . 


TABLE  II 


Comparison  of  PEAPS  and  theoretical  range  increment  values 


D 

Ar  = min (0.57D, 1000) 

Ar  = 5[D] 

50 

28.5 

35.3 

100 

57.0 

50.0 

200 

114.0 

70.7 

400 

228.0 

100.0 

800 

456.0 

141.4 

1600 

912.0 

200.0 

3200 

1000.0 

282.8 

6400 

1000.0 

400.0 

Graphic  outputs  depend  upon  the  Mode  number  chosen. 

Figure  7 is  an  illustration  of  a Mode  0 run  when  CZ  conditions 
are  met.  Figure  8 shows  a Mode  2 run.  These  pictorial 
outputs  may  serve  as  training  aids  in  either  a classroom 
or  a shipboard  environment  to  demonstrate  the  dependence  of 
ray  paths  on  the  vertical  sound  speed  structure.  Mode  1 
has  no  graphic  ray  trace  output. 
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figure  8.  Sample  MODE  2 ray  trace 


IV.  TRANSMISSION  LOSS  GENERATOR 


A.  PATH  SELECTION 

Transmission  loss  is  computed  for  paths  which  meet 
the  following  criteria: 

1.  Surface  Duct  Path  (SDP) 

a.  The  sound  speed  profile  (SSP)  contains  a sonic 

% 

layer.  A sonic  layer  is  defined  as  the  layer  above  the 
first  negative  sound-speed  gradient  (Fig.  9) . 

b.  The  source  and  receiver  are  both  located  within 
the  sonic  layer. 

c.  Sound  remains  trapped  in  the  sonic  layer.  This 
trapping  fails  to  occur  at  very  lov/  frequencies,  where  the 
wavelength  has  become  too  large  to  "fit"  in  the  duct.  Using 
a sound  speed  of  5170  ft/sec  and  a sound  speed  gradient  of 
0.017  sec”1  in  the  sonic  layer,  we  have  from  Urick  [2] 


H 


min 


36.64 


,0.667 
a t 


for  the  minimum  duct  thickness  (H  . ) in  feet  which  will 

nun 

trap  sound  with  wavelength  X in  feet.  Urick  states  that 
this  expression  does  not  represent  a sharp  cutoff.  To 
reduce  programming  steps  for  the  calculator,  the  H 
equation  was  simplified  to 


II  . = 32 

mm 


,0.7 


Table  III  compares  the  frequencies  and  corresponding  minimum 
duct  thickness  as  determined  by  PEAPS , Urick,  and  ASRAP  III 
correction  tables. 

2 . Reliable  Acoustic  Path  (RAP) 

RAP  conditions  always  exist  except  when  the  source 
and  receiver  are  both  located  in  a sonic  layer. 

3 . Bottom  Bounce  Path  (BBP) 

Bottom  Bounce  Path  always  exists. 

4.  Convergence  Zone  Path  (CZP) 

CZP  is  computed  if  an  SSP  with  a sound  speed  excess 
is  present  (Fig.  4) . 

5.  Deep  Sound  Channel  Path  (DSCP) 

The  Deep  Sound  Channel  Path  exists  when  the  source 
and  receiver  are  both  located  in  the  deep  sound  channel  and 
the  sound  speed  profile  can  support  CZ  propagation. 

B.  INPUTS 

The  operator  loads  the  TL/PD  program  and  enters  (from 
the  TL/PD  portion  of  the  worksheet) (Fig.  5)  source  frequency 
(Hz)  , Figure  of  Merit  (dB)  , System  Sigma  (dB)  , wa^e  height 
(ft) , and  MGS  bet tom  type.  Although  the  figure  of  merit 
and  system  sigma  are  used  only  in  the  PD  generator,  they 
are  entered  at  this  time  to  simplify  operator  procedures. 

The  transmission  loss  sub-program  uses  the  following 
values  from  the  ray  trace  generator:  a CZ  indicator,  the 

maximum  and  minimum  CZ  ranges  where  CZ  conditions  exist, 
the  minimum  bottom  bounce  angle,  and  the  depths  and  speeds 
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TABLE  III 


Comparison  of  minimum  duct  thickness 


Frequency 

(Hz)  (ft) 


Minimum  Duct  Thickness  - Hm£n  (ft) 
PEAPS  Urick  AS RAPS  III 


53 

97.5 

789.7 

776.2 

750 

104 

49.7 

492.7 

495.3 

475 

208 

24.9 

303.7 

312.4 

300 

588 

00 

• 

CO 

146.6 

156.2 

150 

1080 

4.8 

95.9 

104.3 

100 

3055 


1.7 


46.4 


52.2 
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of  sound  at  the  critical  points  of  the  sound  speed  profile 
(sonic  layer,  deep  sound  channel  axis,  etc.)- 

C.  TRANSMISSION  LOSS  MODELS 

The  program  looks  at  each  path  selected  and  determines 
transmission  losses  according  to  certain  simple  models. 

Names  used  in  this  section  are  those  used  in  the  FORTRAN 
version  of  the  REAPS  program  (Appendix  D) . 

1.  Surface  Duct  Path  (SDP) 

Figure  10  illustrates  SDP, 

TLSDP  = 201o9(R)  + AR  + DFL*R  for  R < RT  , R < RS 

TLsop  = lOlog(R)  + AR  + DFL*R  + lOlog(RT) 
for  RT  < R < RS 

TLsop  = 10  log (R)  + AR  + DFL*R  + lOlog(RT)  + SSL*R 
for  R > RS  , R > RT 

TLSDP  = 201°9(R)  + AR  + DFL*R  + SSL*R  for  RS  < R < RT 

where  R = the  horizontal  distance  from  the  source  to  the 
receiver  in  yards 

A = the  absorption  coefficient  [3] , 

= + - 40F2  -?]  * 10~3  (dB/yd)  , 

1 + F 4100  + F 

F = frequency  in  kilohertz  (KHz) ; 


DFL  = diffractive  leakage  out  of  the  duct  based  on 
Normal  Mode  Theory  [4  and  5],  given  by 


744f“1,66G~0,33(CD2)“3  * 10“3  (dB/yd),  and 


G 


the  below- layer  gradient  = 


CS3  - CS2 
CD 3 - CD2 


/ 


CD2  = the  sonic  layer  depth  (ft) ; 


CDS  = the  depth  below  CD2  at  which  the  next  positive 
gradient  begins  or  the  bottom  depth  (ft) ; 


CS2,  CS3  = the  sound  speeds  (ft/sec)  which  correspond 
to  CD2  and  CD3,  respectively 

RT  = the  transition  range  (the  range  at  which  the  duct 

is  considered  full  and  cylindrical  spreading  begins) 

m yards  [2]  = 6arccos- (sv/CS2)~ 

SV  = the  speed  of  sound  at  the  source  (ft/sec) . 


RS 


= the  skip  distance  (the  range  between  consecutive 
surface  reflections  of  the  deepest  trapped  ray) 


in  yards 


ro,  _ 2arccos (CS1/CS2) 
1 J 3(R£DINV) 


(yd) 


RADINV  = the  inverse  of  the  radius  of  curvature  of  a 
ray  within  the  surface  duct  = GD/SD. 

SD  = the  depth  of  the  source  of  sound  (ft)  . 

GD  = the  effective  gradient  in  the  sonic  layer 
_ CS2  - CS1 
CD2  - CDl 

0 OB0,71 

SSL  = sea  surface  loss  = ■ — -n-0 (dB/yd) 

Ko 

B = F*H 

H = wave  height  (ft)  . 

The  expression  for  SSL  was  developed  from  a power  regression 
analysis  of  data  presented  in  Reference  6,  and  is  valid  for 
values  of  B between  0 and  30. 
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2.  Reliable  Acoustic  Path  (RAP) 


Front  Fig.  11, 


TL^p  = 201og (K)  + A*K  , 


where  K is  the  slant  range  from  the  source  below  the  layer 
to  the  surface  at  range  R,  in  yards. 

The  maximum  RAP  range  = RL  + RS/2 , where  RL  is  the 
cutoff  range,  defined  as  the  horizontal  distance  from  the 
source  to  the  point  at  which  a ray  just  grazes  the  bottom 
of  the  sonic  layer  (or  surface  if  no  sonic  layer  exists) , 
and  RS  is  the  skip  distance  (described  previously  in  the 
surface  duct  model) . 

3.  Bottom  Bounce  Path  (BBP) 


Four  separate  bottom-bounce  paths  are  possible  between 
the  source  and  receiver:  two  with  one  surface  reflection, 

one  with  two  surface  reflections,  and  one  with  no  surface 
reflections  (Fig.  12) . Incoherent  superposition  of  these 
rays  at  the  receiver  results  in  a 6 dB  reduction  in  the 
transmission  loss.  For  emission  angles  greater  than  about 
15°,  the  radius  of  curvature  of  a ray  is  usually  large  enough 
to  assume  straight-line  propagation  of  the  rays  from  the 
source  to  the  bottom  and  back  to  the  surface.  A ray  from  the 
source  will  strike  the  bottom  at  an  angle  of  incidence 
estimated  by 
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THBL  = tan 


rCD6  - SD, 
1 XB  J 


where 


XB 


R (CD6  - SD) 

N T2Cl5'6  - SD) 


(yds) 


is  the  horizontal  distance  between  the  source  and  the  range 
of  the  first  bottom  bounce,  CD6  is  the  bottom  depth  (ft) , 
and  N is  the  number  of  bottom  reflections. 

The  distance  traveled  by  the  four  rays  which  strike 
the  bottom  once  can  all  be  approximated  by 


cos (THBL)  (yds) 

for  reasonably  small  source  and  receiver  depths  or  large 
range.  Therefore,  an  approximate  transmission  loss  can  be 
written  as 

P A*R 

TLBB  = 201o9  ^cqs (THBL) ^ + cos (THBL)  + N*BL  “ 6 

where  the  subtraction  of  6 takes  into  account  the  6 dB  gain 
caused  by  multiple  bottom  bounce  paths,  and  BL  is  a bottom 
loss  term  which- was  determined  through  regression  analysis 
of  bottom  loss  data  presented  in  Reference  7.  Specific  values 
for  BL  are  determined  by  a sorting  routine  which  uses  the 
input  frequency  and  bottom  type.  THBL  is  the  angle  at  which 
the  ray  strikes  the  bottom. 


40 


4 . Convergence  Zone  Path  (CZP) 


See  Fig.  4. 

TLCZ  = 201og (R)  + A*R 

when  R is  outside  the  CZ  annulus  and 


TLCZ  = 201og (R)  + A*R  - GCZ 

when  R is  bounded  by  CZSTRT  and  CZSTOP . The  CZ  annulus  is 
defined  as  the  width  of  the  zone,  CZSTOP  - CZSTRT,  in  yards. 
GCZ  is  the  reduction  in  transmission  loss  caused  by  the 
increase  in  sound  pressure  level  (SPL)  which  results  from  ray 
focusing  within  the  ray  family  forming  the  CZ  annulus. 


RANGE  = R - CZM1N  (yds) 

X = CZ  annular  width  (yds) 

GCZ MAX  = 15.4545  - 4.545  * 10"4X  (yds) 


The  equations  for  GCZ  and  GCZMAX  are  empirical  and  were 
derived  from  analysis  of  numerous  transmission  loss  curves 
obtained  as  output  from  the  P'ACT  model  [5]  . The  analysis 
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established  that  the  gain  distribution  within  a convergence 
zone  closely  resembled  a Rayleigh  distribution  and  that  the 
gain  tended  to  a maximum  at  approximately  18%  of  the  annular 
width.  Thus,  GCZ  was  modelled  as  a Rayleigh  distribution 
with  the  peak  occurring  at  a distance  from  CZSTRT  equal  to 
18%  of  the  annular  width  (Fig.  13) . 

5.  Deep  Sound  Channel  Path  (DSCP) 

See  Fig.  14. 

TLDSCP  = 101o9(R)  + A*R  + lOlog(RTP) 
where  RTP  is  the  deep  sound  channel  transition  range, 

RTP  -SWH  ,§,* 

D = CD3  - CD2  in  feet, 

CD3  = depth  of  the  sound  channel  axis  (ft) , 

Z = SD  - CD2  (ft) , 

CZMIN  is  in  yards. 

The  transmission  loss  term  is  valid  only  for  R > RTP, 
although  it  is  used  for  all  R when  this  path  is  present. 

D.  COMBINED  TRANSMISSION  LOSS 

Since  transmission  loss  caused  by  bottom  bounce  is  always 
present,  the  total  transmission  loss  (STL)  is  initialized 
with  TLbb  after  each  increment  of  range.  Transmission  loss 
for  subsequent  propagation  paths  meeting  their  respective 
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figure  13.  Convergence  Zone  Gain  Model 
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existence  criteria  are  computed  for  the  same  range  and 
combined  with  STL  as  follows:  If  © signifies  the  combining 

of  decibels  (dB)  for  independent  (phase  incoherent)  signals, 
then 

STL  © TL (i)  = STL  - A 

where  TL(i)  is  the  transmission  loss  for  the  specific 
propagation  path  (CZ,  SDP,  DSC,  RAP,  as  they  exist). 

A = lOlog (1 + 10“°*1J) 

J = jsTL  - TL  (i)  . 

E.  OUTPUT 

The  program  plots  the  combined  transmission  loss  (STL) 
in  dB  versus  range  in  nautical  miles  (NM)  , with  a reference 
line  for  the  input  figure  of  merit  (Fig.  15) . 

F.  SHORTCOMINGS  OF  THE  MODELS 

Since  the  models  were  developed  for  use  in  a relatively 
small  programmable  desk-top  calculator,  limitations  in  the 
available  storage  and  the  number  of  program  steps  have 
resulted  in  certain  shortcomings: 

1.  Phase  coherence  is  neglected. 

2.  Sea  surface  scattering  is  not  taken  into  account  in 
calculating  TL  for  the  convergence  zone  model. 

3.  If  the  source  is  in  the  sonic  layer  and  the  duct 

thickness  is  less  than  II  . , propagation  in  the  layer  is 

nun 

neglected. 
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4.  The  use  of  straight-line  propagation  will  induce 
some  error  in  the  angles  computed  for  bottom  bounce 
parameters  and  in  the  transmission  loss  computation. 

5.  Convergence  zones  are  considered  when  any  sound 
speed  excess  is  present. 

6.  The  convergence  zone  gain  term,  GCZ,  is  found  from 

an  empirical  formula:  . •>.re  is  no  physical  justification 

for  the  use  of  the  Rayleigh  distribution. 

7.  The  model  for  the  deep  sound-channel  overestimates 
transmission  loss  for  ranges  less  than  the  (deep  sound- 
channel)  transition  range . 

8.  For  reliable  acoustic  path,  steep-gradient  curvature 
is  not  considered;  transmission  loss  is  computed  for  straight 
line  distances  between  the  source  and  receiver. 

9.  The  absorption  loss  term  uses  Thorp's  formula  for 
all  frequencies.  As  shown  in  Fig.  16,  the  formula  is  good 
for  frequencies  greater  than  20  Hertz. 

10.  The  diffraction  loss  term,  DFL,  is  the  same  term  used 

in  FACT  [5].  The  following  is  a quote  from  Reference  5: 

A preliminary  analysis  of  this  model  (Labianca,  1972) 
has  shown  that  the  leakage  term  does  not  correspond  to 
the  mechanism  which  dominates  duct  leakage  and,  in 
particular,  has  the  wrong  dependence  on  the  below-layer 
(thermocline)  gradient.  An  intensive  investigation  of 
all  available  surface  duct  models  has  been  completed 
recently  by  NUC  (Naval  Undersea  Center)  for  LRAPP  (Long 
Range  Acoustic  Propagation  Project) (Morris,  1974) . This 
study  compared  several  models  with  each  other  and  the 
SUDS  (Surface  Duct  Sonar)  transmission  loss  data 
(Cummins,  1972)  . Since  none  of  the  candidate  models 
emerged  as  clearly  superior,  AESD  (Acoustic  Environmental 
Support  Detachment)  is  currently  studying  possible 
modifications  to  the  surface-duct  equations.  [5] 
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FREQUENCY  IN  KHZ 


Measured  data  taken  from  the  North  Atlantic  (Thorn  and 
3ernier,  1959  and  1962;  Susswan,  MacDonald,  and  Kanahis, 

1963) ,  Mediterranean  (i.al  lenient  and  Waterman,  1963;  leroy 

1964) ,  and  Northern  Pacific  (Greer  and  Polam,  1954; 

Sheehy  and  Halley,  1957)  [3j. 


Thorp's  Ea.:-  Attenuation 


0.  IF' 


H-40F' 


x 10‘J  dB / 


Figure  16 


11.  Convergence  zone  parameters  are  determined  at  the 
surface.  Errors  exist  when  source  and/or  receiver  are  at 
appreciable  depths. 

12.  When  CZ  conditions  are  not  present,  the  model 
assumes  there  is  no  deep  sound  channel. 

13.  Leakage  of  sound  into  a shadow  zone  is  not  con- 
sidered. If  there  are  no  CZ  conditions,  then  the  only 
illumination  of  the  shadow  zone  comes  from  bottom  bounce. 

If  CZ  conditions  exist  the  shadow  zone  is  illuminated  by 

a path  having  spherical  spreading  in  addition  to  the  bottom 
bounce  illumination. 
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V.  PROBABILITY  OF  DETECTION  (PD)  GENERATOR 


A . INPUTS 

The  sub-program  for  probability  of  detection  contains 
values  from  the  cumulative  distribution  function  (CDF)  for 
a standardized  normal  random  variable,  CDF  = F(x)  , Values 
of  x corresponding  to  F(x)  = 0.5,  0.55,  ...,  0.95  are  stored. 

The  figure  of  merit  (FOM)  and  system  sigma  are  entered 
by  the  operator  as  part  of  the  transmission  loss  initializa- 
tion procedure.  System  sigma  is  defined  as  the  square  root 
of  the  sum  of  the  variances  of  the  individual  terms  in  the 
passive  sonar  equation: 

System  sigma  = (a|L  + + c ^ + o^)h 

where  o is  the  deviation. 

Rule-of-thumb  values  for  system  sigma  are  included  in 
Operator  Note  3 (Appendix  C)  for  those  cases  in  which  variances 
of  the  individual  term  are  unknown. 

B.  MODEL 

The  detection  variable  X is  defined  as  the  ratio  of  the 
signal  excess, 

SE  = FOM  - TL  , 

to  the  system  sigma: 

x = SE 

System  Sigma 


58 


The  computed  X is  compared  to  the  stored  values  of  x 
to  determine  probability  of  detection  at  which  PD  = F(X) . 

C.  OUTPUT 

The  program  plots  a probability  of  detection  (PD)  versus 
range  directly  below  the  TL  versus  range  plot  (Fig.  15) . 

The  PD  routine  computes  PD  in  increments  of  0.05  so  that  a 
smooth  curve  does  not  result. 


VI . CONCLUDING  REMARKS 


Advanced  development  of  an  on-board  environmental 
prediction  system  resulted  in  the  Integrated  Command  ASW 
Prediction  System  (ICAPS)  computer  program.  ICAPS  uses  the 
FACT  transmission  loss  program  so  that  it  requires  computer 
resources  comparable  to  those  available  on  aircraft  carriers. 
This  development  effort  is  continuing  with  emphasis  upon  the 
development  of  a smaller  version  of  FACT.  A viable  alternative 
would  be  an  augmented  version  of  PEAPS. 

Fleet  introduction  of  an  environmental  prediction  system 
need  not  be  delayed  until  a foolproof  system  is  available. 

With  PEAPS  available,  even  in  its  present  limited  form, 
operational  evaluation  of  an  interim  system  with  some  known 
shortcomings  is  entirely  feasible.  A continuing  dialogue, 
carried  on  between  the  developmental  agency  and  the  operating 
forces,  would  facilitate  implementation  of  system  improvements 
based  upon  fleet  experience. 

The  operating  forces  have  never  had  a tactical  environ- 
mental prediction  capability  available  at  the  individual 
platform  level.  PEAPS  provides  an  inexpensive  way  to  evaluate 
the  operational  utility  of  such  a system. 
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APPENDIX  A 


PEAPS  HARDWARE  DESCRIPTION 

For  this  program  the  Hewlett-Packard  Model  9810A 
calculator  must  be  configured  with  the  following  options: 
Option  001:  Total  of  111  data-storage  registers; 

Option  003:  Total  of  2036  program  steps; 

Option  004:  Printer. 

Read-only-memory  units  (ROM's)  required  are: 
Mathematics  ROM  - hp  - 112 10 A 
Printer  Alpha  ROM  - hp  - 11211A 
Plotter  ROM  - hp  - 11215A. 

The  Hewlett-Packard  Model  9862A  Plotter  is  required  to 
produce  the  graphical  output. 


61 


PEAPS  OPERATOR  ’WORKSHEET 


(See  PEAPS  Operator  Note  Sheets) 
BT  Information 


DTG: 


Location: 


1.  Turn  the  calculator  and  plotter  on.  "LOAD"  the  SSP/Ray 
trace  program  cards.  Put  the  plotting  paper  on  the 
plotter  and  set  the  "GRAPH  LIMITS"  of  the  plotter. 

2.  Press  "STOP,  END,  CONTINUE"  (this  sequence  is  used  when 
a new  SSP  is  to  be  generated. 

Note  1;  If  you  enter  an  erroneous  value: 

a.  If  you  have  not  pressed  "CONTINUE",  press 
"CLR  X",  then  enter  the  correct  value. 

b.  If  you  have  pressed  "CONTINUE",  you  must 
re-enter  all  the  points  again;  press  "GO  TO, 
LABEL,  A.",  then  continue  with  Step  3 below. 

3.  Press  "CONTINUE". 

4.  SSP  inputs:  ENTER  number  of  BT  points  (15  maximum) 

Press  "CONTINUE"  after  each  value  is  entered  to  advance 

the  program  Output 

Input  Input  Input  Speed  of  Sound 

Point  Depth  (ft)  Temp  (°?)  Salinity  (ppt)  (ft/sec) 


Note  2 : Change  the  plotting  paper  on  the  plotter.  Enter 

the  ray  trace  program  MODE  indicator. 

0 = Automatic  compute  and  plot 

MODE:  1 = Automatic  compute,  no  plot 

2 = Manual  inputs  for  plotting.  Automatic 
compute . 

If  a ray  trace  input  error  is  made,  press  "GOTO, 

LABEL,  B,  CONTINUE",  then  input  the  MODEL  indicator 
and  ray  trace  data  again.  (This  sequence  is  used 
to  generate  a new  trace  using  the  old  SSP  data) 

Press  "CONTINUE" 

\ 

4.  Ray  trace  inputs: 

Press  "CONTINUE"  after  each  input  to  advance  the  program. 

Source  Depth:  (ft) 

If  Mode  2 is  used: 

Minimum  Ray  Angle  x 2:  (degrees) 

Angle  Increment  x 2:  (degrees) 

Maximum  Ray  Angle  x 2: (degrees) 

Note  3:  If,  while  in  Mode  0 and  after  the  bottom  limiting 

ray  ha  - bee,;  drawn,  a channelled  ray  is  drawn  to 
the  plot's  range  limit;  see  Operator  Note  sheet 
section  Note  2E. 

Change  the  plotting  ..per  on  the  plotter.  Press 
"STOP,  END".  LOAD  tns  TLOSS/PD  program  cards. 

Press  "STOP /END" . (The  "STOP,  END"  sequence  is 
used  when  a TLOSS/PD  input  error  is  made,  or  a new 
TLOSS/PD  plot  is  to  be  generated  based  on  the  old 
SSP/Ray  trace  data.) 

5 . TLOSS/PD  inputs : 

Press  "CONTINUE"  after  each  input  to  advance  the  program. 
Press  "CONTINUE" 

Source  Frequency:  (Hz) 

Figure  of  Merit  (FOM)  : (aB) 

System  Sigma:  (use  10  if  unknown) 

Wave  Height:  #123456789  10  11  12  (ft) 

Receiver  Depth:  (ft) 

Bottom  Type:  123456789 
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If  a new  ray  trace  is  to  be  generated  using  the 
old  SVP  data  after  the  TLOSS/PD  program  has  been 
BH,  Press  "STOPTeND"  . LOAD  the  SVP/Ray  trace 

program  cards. 

Press  "GOTO,  LABEL,  B,  CONTINUE",  then  see  Note  2 
above . 


is  inserted  into  the  calculator  as  follows  (it  is  assumed 
the  calculator  is  turned  on  and  paper  has  been  placed  on 
the  plotter) . 

1.  Press;  END,  LOAD. 

2.  Insert  the  cards  starting  with  side  1 sequentially 
through  side  4. 

3.  Press;  END,  CONTINUE. 

After  the  SSP  axes  have  been  drawn  the  printer  will  indicate 
"SEE  NOTE  1*',  referring  the  operator  ot  Note  1. 

Note  1 

A.  The  sound  speed  profile  (SSP)  generator  will  accept 
up  to  fifteen  (15)  points.  These  points  are  entered  from 
the  PEAPS  Computer  Operation  Worksheet  starting  with  the 
first  point  being  the  surface  (depth  = 0)  temperature  and 
| salinity.  The  data  are  entered  as  follows i 

1.  Press;  CONTINUE 

The  printer  will  request  "NO.  PTS".  Enter  the 
number  of  BT  points  to  bo  entered  as  indicated  on  the 


worksheet. 


2.  Press:  CONTINUE . 


[ 

r 


P 
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The  printer  will  request  DEPTH.  Enter  the  depth 
in  feet  of  the  point  for  which  data  are  being  entered  (enter 
0 for  point  1) . 

3.  Press:  CONTINUE. 

The  printer  will  print  the  depth  entered  and 
request  TEMP.  Enter  the  temperature  in  °F  for  the  point 
of  interest. 

4.  Press:  CONTINUE. 

The  printer  will  print  the  temperature  entered 
and  request  SALIN.  Enter  the  salinity  in  parts  per  thousand 
for  the  point  of  interest. 

5.  Press:  CONTINUE. 

The  printer  will  indicate  the  salinity  entered 
and  the  speed  of  sound  in  feet  per  second  corresponding  to 
the  depth,  temperature,  and  salinity  values  immediately 
preceding  the  sound  speed  value.  The  plotter  will  plot  the 
speed  of  sound  indicated  on  the  printer  at  the  corresponding 
depth.  If  the  data  entered  are  from  a point  less  than  the 
total  number  of  points  to  be  entered,  the  printer  will  indi- 
cate DEPTH.  Enter  the  depth  of  .he  next  point  of  interest 
and  continue  from  step  3 above.  If  the  date  entered  is  for 
the  last  point,  the  program  will  discontinue  plotting  and 
the  printer  will  indicate  "SEE  NOTE  2".  This  refers  to 
Operator  Note  2 below. 

B.  If  the  BT  trace  does  not  reach  the  bottom: 
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1.  Use  at  least  two  points  from  historical  data 


and  enter  these  as  the  last  two  points.  The  more  points 
entered,  the  more  accurate  the  plot  if  the  bottom  is  con- 
siderably deeper  than  the  BT  trace  depth.  See  Section  G 
below  for  merging  data. 

2.  If  historical  data  are  not  available,  enter  the 
last  BT  temperature  as  the  bottom  temperature,  the  bottom 
depth  as  measured  or  approximated,  and  35  ppt  as  salinity. 

C.  If  salinity  values  are  not  available  from  measurements 
or  historical  data,  enter  35  ppt. 

D.  If  you  enter  an  erroneous  value: 

1.  If  you  have  not  pressed  CONTINUE,  press  "CLR  X"  , 
then  enter  the  correct  value. 

2.  If  you  have  pressed  CONTINUE,  you  must  reenter 
all  the  points  again;  press  "GO  TO,  LABEL,  A",  then  continue 
from  step  2 in  Section  A above. 

E.  If  you  have  finished  with  one  sound  speed  profile 
and  want  to  enter  a new  one,  press  "END",  then  go  to  step  1 
in  Section  A above. 

F.  Pressing  "CONTINUE"  after  each  input  advances  the 
program  and  enters  the  input  point  into  storage. 

G.  If  the  .in-situ  and  historical  data  do  not  match 
closely  together,  merging  becomes  necessary.  There  is  no 
single  correct  method  for  merging  historical  and  in-situ 
data.  Suggested  procedures  arc  to  fare  (draw  a smooth  curve) 
the  two  profiles  together,  to  extrapolate  the  in-situ  trace 
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down  to  the  historical  data,  or  to  extend  the  profiles  in 
the  area  of  interest  until  they  meet.  Fig.  17a  shows  a plot 
of  in-situ  SSP  and  historical  data  on  the  same  axis.  It  is 
suggested  that  this  be  done  for  all  SSP  plots  if  historical 
data  are  available.  A trace  of  the  historical  data  need 
not  start  at  the  surface.  A depth  500  feet  above  the  last 
BT  depth  should  suffice.  Figures  17b  and  17c  show  examples 
of  extending  the  two  traces  to  achieve  a merge.  In  any 
event,  a smooth  blending  of  the  traces  should  be  attempted. 
In  some  instances  it  may  be  necessary  to  try  more  than  one 
technique  in  order  to  achieve  a reasonably  correct  facsimile 
of  the  sound  speed  structure. 

Note  2 

This  Note  deals  with  procedures  for  the  ray  trace 
generator . 

A.  Change  the  plotting  paper.  , 

B.  Enter  the  ray  trace  program  Mode  indicator  (0,  1, 
or  2)  . 

Mode  0 gives  an  automatic  search  for  and  plotting 
of  convergence  zone  critical  rays.  The  minimum  source 
emission  angle  that  will  produce  a bottom  bounce,  and  the 
minimum  (inside)  and  maximum  (outside)  convergence  zone 
annulus  ranges  (if  CZ  is  present)  are  stored. 

Mode  1 gives  the  same  computation  as  Mode  0 except 
that  the  ray  traces  are  not  plotted. 
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Mode  2:  the  operator  enters  twice  the  minimum  and 

maximum  emission  angles  to  be  traced  from  the  source  along 
with  twice  the  incremental  angle  between.  Doubling  the 
desired  angles  and  angle  increment  is  necessary  since  the 
program  internally  halves  all  input  angles  prior  to  conver- 
sion to  radians  and  storing.  When  the  manual  input  plotting 
is  completed,  the  program  changes  to  Mode  1 to  determine 
the  critical  parameters  to  be  used  by  follow-on  programs. 

C.  Program  Mode  0 is  intended  to  be  the  primary  Mode 
for  fleet  use.  Mode  1 is  intended  for  use  when  a ray  trace 
is  not  desired.  Mode  2 is  intended  to  be  used  for  ray  trace 
plots  of  specific  source  emission  angles  as  well  as  for 
determination  of  critical  parameters  for  follow-on  programs  - 
it  is  foreseen  as  a training  aid. 

D.  If  you  have  made  an  incorrect  input  and  entered  it 

with  the  CONTINUE,  press:  GO  TO,  LABEL,  B,  CONTINUE. 

Re-enter  the  ray  trace  input  parameters  starting  with  the 
Mode  indicator. 

E.  It  is  possible,  under  certain  sound  velocity  profile 
conditions,  for  the  ray  trace  program  in  the  automatic  modes 
to  trace  a ducted  ray  to  the  end  of  the  plot  while  searching 
for  the  convergence  zone  critical  ranges.  If  this  situation 
should  occur  you  must  check  the  stored  convergence  zone 
ranges  as  follows: 

1.  Press:  X<-(  ),  008,  1 , 2000,  t. 
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This  calls  out  the  minimum  convergence  zone 


range  and  converts  it  from  yards  to  nautical  miles. 

2.  Compare  the  value  in  the  Y-register  with  the 
minimum  convergence  zone  range  indicated  on  the  ray  trace 
plot.  If  they  differ  considerably, 

a.  Enter  the  minimum  convergence  zone  range 
from  the  ray  trace  plot; 

b.  Press:  t-,  2000,  x (multiply)  , Y-»-(  ) , 008. 

This  converts  the  ray  trace  plot  minimum 

CZ  range  from  nautical  miles  to  yards  and  stores  it  in  the 
appropriate  storage  register. 

3.  Press  X«-(  ),  009,  \ , 2000,  *. 

This  calls  out  the  maximum  CZ  range  and  converts 
it  to  nautical  miles. 

4.  Compare  the  value  in  the  Y-register  with  the 
maximum  CZ  range  off  the  ray  trace  plot.  If  they  differ 
considerably; 

a.  Enter  the  maximum  CZ  range  from  the  ray 

trace  plot. 

b.  Press:  f,  2000,  X (multiply),  Y-*-(  ),  009. 

This  converts  the  ray  trace  plot  maximum 

CZ  range  from  nautical  miles  to  yards  and  stores  it  in  the 
appropriate  storage  register. 

5.  Proceed  with  the  transmission  loss/probability 
of  detection  program. 
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"TH|  |T*» 


F.  If  you  want  to  do  another  ray  trace  and  plot  using 
the  same  sound  speed  profile  but  different  input  parameters: 
1.  Press:  GO  TO,  LABEL,  B,  CONTINUE. 


2.  Enter  the  new  ray  trace  inputs  starting  with  the 
Mode  indicator. 

G.  ''SEE  NOTE  3"  and  "END"  will  be  printed  on  the  output 
tape  and  the  calculator's  X,  Y , and  Z registers  will  illuminate 
when  the  ray  trace  portion  of  the  program  is  completed. 

H.  "NO  CZ",  "SEE  NOTE  3",  and  "END"  will  be  printed  on 
the  output  tape  to  indicate  that  the  program  has  ended  and 
there  is  no  convergence  zone  present. 

I.  If  "NO  CZ" , "SEE  NOTE  3",  and  "END"  is  printed  out 
and  you  would  like  a ray  trace  plot,  you  will  have  to  use 
Mode  2 and  manual  emission  angle  parameters  to  get  a plot 
(see  E above) . At  the  end  of  the  Mode  2 plots  "NO  CZ", 

"SEE  NOTE  3",  and  "END"  will  be  printed  on  the  output  tape. 

J.  If  another  sound  velocity  profile  is  to  be  entered 

for  subsequent  ray  trace  plots,  press:  END,  CONTINUE,  to 

start  the  entire  sequence  again. 


TRANSMISSION  LOSS/PROBABILITY  OF  DETECTION  PROGRAM 
Note  3 

A.  Change  the  plotting  paper. 

B.  The  TL/PD  magnetic  cards  are  loaded  in  the  same 
manner  as  the  SSP/Ray  Trace  cards. 

C.  Press:  STOP,  END  to  initialize  the  program  function. 

This  sequence  is  used  when  a TLOSS/PD  entering  error  is  made, 
or  a new  TLOSS/PD  plot  is  to  be  generated  based  on  the  same 
SSP/Ray  Trace  generated  data. 

D.  Source  frequency  is  to  be  entered  in  Hertz  (Hz) . 

E.  Figure  of  Merit:  ensure  that  you  use  all  the  same 

units  of  reference  (pPa,  or  pb)  . DO  NOT  add  100  to  the 
computed  FOM. 

F.  System  Sigma 

If  the  variances  of  the  sonar  equation  terms  are 
known,  enter  the  square  root  of  the  sum  of  the  variances. 

If  the  variances  are  unknown,  choose  one  of  the  following: 

1.  Use  6 dB  if  ambient  noise  measurements  have  been 
made,  and  submarine  speed  and  type  are  known. 

2.  Use  0 dB  if  ambient  noise  is  estimated  from  a 
forecast,  submarine  speed  is  known  to  within  3 kts,  and 
submarine  type  is  known . 

3.  Use  10  dB  if  ambient  noise  is  estimated  from 
forecasts  and  submarine  speed  and  type  are  uncertain. 
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G.  Wave  Height  (IWH) 


Sea  State 

Term 

Height  of  Wave  (ft) 

1 

Calm 

0 to  1 

2 

Wavelets 

1 to  2 

3 

Slight 

2 to  4 

4 

Moderate 

4 to  8 

5 

Rough 

8 to  13 

6 

Very  rough 

13  to  20 

7 

High 

20  to  30 

8 

Very  high 

30  to  45 

5 

Phenomenal 

over  45 

H.  A bottom  type 

indicator  (MGS)  of 

1 indicates  little 

sound  loss;  a bottom  type  indicator  of  9 

indicates  a large 

sound  loss. 

I.  If  new  TLOSS/PD  data  is  to  be  entered  to  make 
additional  plots  from  the  same  ray  trace  data,  press  STOP, 
END.  Change  the  plotting  paper,  then  press  CONTINUE.  Then 
enter  the  new  data. 

J.  If  a new  ray  trace  is  to  be  generated  from  the  old 

SSP  information  after  the  transmission  loss/probability  of 
detection  program  has  been  run:  press:  STOP,  END,  LOAD. 

Load  the  SSP/Ray  Trace  program  using  the  magnetic  cards, 
change  the  plotting  paper,  press:  GO  TO,  LABEL,  B,  CONTINUE, 

then  follow  the  instructions  in  Note  2. 
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APPENDIX  D 


Flow  charts,  program  listings,  and  output  from  the 
HP9810A  and  IBM  360  are  included  within  this  Appendix. 


Speed  of  Sound  Profile  / Ray 
Trace  Program 


LABEL  A 
See  Note  1 


Input  the  # 
of  SSP  pts* 


■w 


(2) 


Input  Depth 
(ft) 


Input  (°P) 
Temperature 


Input  ( °/ oo  )j 
Salinity 


Output  Sound 
Speed  & plot 
vs  depth 


z 

Draw  the 

Sort  out  and 

Speed  of 

store  the 

Sound 

critical  Speed 

Profile  

of  Sound 

k^Axe^|-^  (1) 

Profile  points 

(3) 


, have 
all  the  SSP 
joints  beetle*, 
nput'? 


a. 


W 


jABEL  b 
See  Note  2 


Input  MODE 
indicator 


Initilize. 
the  needed 
parameters 
for  Ray 
Tracing 


Input  Source) 
Depth  (ft) 


MODES i 

0 
1 
2 


Auto  compute  and  plot 
Auto  compute,  No  plot 
Manual  input  of  ray 
family  for  plotting, 
then  MODE  i Auto 
compute  with  no  plot 
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(2)  Sample  Speed  of  Sour.d 
Profile  program  input  / 
printer  output 


IfF.P  I'H 

0 . u 0 * 

DEPTH 

400 , 00* 

'DEPTH 

1 1 00 .110* 

••  EI-IP 

5 8 . 6 0 * 

TEMP 

50 . 50* 

TEMP 

45.30* 

SrlLIH 

35 . oo* 

•SOL  111 

35. 00* 

SRLIH 

35 . UO* 

SPii 

4341 . ’■  * 0 

S p II 

4838 . 36 

SPD 

4877.81 

1'h.r'  i H 

i 20. 110* 

:epth 

650 . uO* 

DEPTH 

1 300 . 00* 

"EilP 

53. 1 0* 

TEMP 

47 . 30* 

’PIP 

45.00* 

'.■Hi.  I H 

35  00* 

SflLHI 

3 5 . 0 0 * 

8HL1H 

35 . 00* 

4343 . 3-4 

s pp 

4835. 12 

SPD 

4875.51 

•TH 

175. UO* 

DEPTH' 

870 . 00* 

DEPTH 

1380.OO* 

i.HP 

57  UO* 

TEMP 

47.20* 

VL'MF 

44  _ 50* 

SftLJ  H 

35 . UO* 

8ft L. ' M 

33 „ OO* 

• : ML  IN 

35 . O0* 

ti 

4335. • 5 

3 PH 

4884 . 56 

SPD 

4873.26 

DEPTH 

liEF'  1 H 

DEPTH 

200  OO* 

300 . 00* 

1 500 . 08* 

7 EMI 

TEMP 

TMP 

56 . UO* 

47.70* 

43.70* 

SLLIH 

SHLIH 

Sf’iL  1 n 

35 . 00* 

35 . 00* 

35 . 00* 

SF  > i 

SI  D 

SPD 

4329.50 

4887  "t. 

4863 . *»9 

.TPT  H 

DEPTH 

DEPTH 

225.U0* 

1 075 . UO* 

7500 . uo* 

TEMP 

TEMP 

‘TTIP 

54 . HO* 

46. 50* 

43. "0* 

•-.ML  1 H 

2 Mi.  J n 

siTi  nr 

35 . H0* 

35 . UO* 

35 . uO* 

•MU 

SPD 

ypD 

4917. '"8 

4 UK  3 ' ' t* 

4 9 6 . U 

see:  mote  2 
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(4)  SSP  critical  points 
selection  program 


(4)  SSP  critical  points  selection  program 


(4)  SSP  critical  points  selection  program 


Speed  of  Sound  Profile  / Ray  Trace 
Program 


RRNGE  CNM3 


B j Speed  of  Sound  Profile  / Ray  Trace 
Program 


MODE  = 2 

v ? . 


/THMIN=THPP\ 
/set  plotting' 
i indicator  , 


Trace  the 
Ray 


the  ray  a 
r 
? 


THMIN=THMIN- . 1 ' 

clear  plotting 
indicator 


set  FLAG 


rRBB-TRMIN\ 
set  ICfc  > 
set  plotting 
indicator  / 
CZmax-R  / 


the  FLAG 


Draw  the  ray 
that  just 
graz.es  the 
bottom  


THP»  se 
lotting 
indica- 


THMIK=THMIN+. 5 
set  plotting 
indicator 


(8)  Ray  Trace  Program 


Call  Point 


(8)  Ray  Trace  Program 
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RFtNEB:  CNKO 


CWKJ3  253N 


RFINBS:  CNM3 


RFJiSKE  CNM 


CNM3 


0000-- 
0001  — 
0002"" 
0003“* 
0004“' 

0005- - 

0006- 

0007- 

0008- 

0009- 

0010- 
0011- 
0012- 
0013“ 

0014- 

0015- 

0016- 

0017- 

0018- 

0019- 

0020- 
0021- 
0022- 

0023- 

0024- 

0025- 

0026- 

0027- 

0028- 
0029- 
0030* 
003 1 ■ 
0032* 
0033- 
0034' 

0035 

0036 

0037 

0038 
'0039 
0040 
004  1 


004 


GTO — 

- 4 4 

0O54  — 

4 — 

-04 

O055-- 

F -- 

-85 

0O56-- 

YTO— 

-40 

8057— 

I HD  — 

-31 

0058-- 

b -- 

-14 

0O59-- 



-03 

0060-“ 

XT  0 — 

-23 

0061 

-34 

0062— 

o.  — 

-13 

0863-- 

6 -- 

-06 

0 064-- 

XTO  — 

-23 

0065-- 

— 

-34 

O O 66 — 

b -- 

-14 

0867-- 

-::fr— 

— l 

8068-- 

-I  HD— 

-31 

OO 69-- 

- o - - 

-13 

087O- 

- X TO-- 

-23 

0071- 

-I  MB  — 

-31 

8072- 

- b -- 

— 14 

8073- 

- 3 -- 

—03 

0074- 

-XTQ  — 

" C-’-' 

0075- 

- + -- 

—33 

0876- 

- oi  —13 

0077- 

"•t 

- f -■ 

—87 

8079- 

-XTO— 23 

8079- 

- + —33 

0088- 

- b - 

— 14 

0081- 

- 1 - 

— 01 

8082- 

-XTO- 

--23 

8083- 

- + ~ 

—33 

O084- 

—82 

0085- 

-S/R- 

— » 

t » 

0886- 

-XFR- 

_ _ "7 

8887- 

-I  HD- 

—3 1 

0888- 

— — 

— 13 

0889- 

- UP- 

-5-V 

w » 

0898- 

--02 

009  1 - 

-XTO- 

— 23 

0892- 

•!-  -• 

— 3 y 

8893- 

— — * 

— 13 

0094- 

—XFR- 

--67 

0895- 

— I MB- 

—31 

0096- 

- - ft  ~ 

— 13 

0897- 

0044  — 8 ■'  R 77 

FMT— 42 


0045 — 
0 0 46  — 

0047— 

0048- - 
•0049 — 
0058  — 
005 1 — 


1 —01 
8 10 

0  OO 

5 05 

2 02 

1 —01 

0 ---00 
UP 27 


8098- 

0099- 

0100- 
0 101' 
0102 

0103 

0104 

0105 

0106 
0107 


4  04 

A 06 

0 00 

0 00 

FMT 42 

1 —01 
2 -—02 
2 02 

5  05 

0 00 

0 00 

1JP 27 

CHS 32 


EEX — 

-26 

4 — 

-84 

FMT  — 

-42 

1 -- 

-01 

o 

V 

-03 

0 

-00 

UP  — 

-27 

4 __ 

-04 

“7  _ _ 

1 

-07 

0 

-00 

- 0 -- 

-00 

-FMT  — 

-42 

- 1 — 

-01 

- UP  — 

-27 

. ! — 

-01 

- 0 

-00 

-FMT  — 

--42 

- 1 ---01 

- 5 05 

- 5 85 

- 0 00 

— n — 00 

-FMT—- -42 


l 

— 01 

6 

06 

• 

21 

1 

— 01 

0 

08 

FMT 

---42 

1 

01 

1 

0] 

1 

01 

r> 

w 

02 

1 

— 87 

5 

—05 

UF 

— — — ■ fci  r 

4 

—04 

- 9 

— ] 1 

- 0 

00 

- 0 

00 

0108—  FMT— 42 

0109—  1 01 

0110—  UP  — 27 

0111—  FMT— 42 

0112—  1 —01 

01 13— FMT 42 

01 14  — V T 0 - t — 4 0 

0115—  A' St- 

Oll  6--  E —60 
8117—  E —60 

0118 — B 63 

81 19— CUT 47 

0128— K<Y 52 

0121 — F 1 6 

8122— XTO  — 23 

0123—  HIV 35 

0124- -YT0 48 

81 25--  E 60 

01 26--  C 61 

8127  — CUT  — 47 
0128—  YE— 24 
8129— CUT ‘47 


0130  — 

1 - 

— 01 

0131  — 

0 * 

00 

0132-- 

0 ' 

00 

0133 — P 

■SE* 

— 57 

0134  — FMT- 

42 

8135— 

4 ■ 

04 

0136 — 

6 

06 

0137  — 

9 

11 

0138  — 

0 

80 

0139 — t 

iTO 

— — — 23 

O140  — 

ft 

1 3 

0141  — 

4 

04 

O 1 42-- 

1* 

07 

8143--: 

<T0 

23 

0144  — 

b 

14 

0145  — 

3 

03 

0146  — 

0 

08 

814?  — 

0 

0 0 

0148— 

UP 

—7 

1 “ •“  c. 

81 49 — 

ft 

13 

015O— FMT 42 

8151—  1 —01 

81 52- -  UP 27 

0153—  b 14 

0154—  FMT— 42 

0155—  1 01 

8 1 56— PUT 45 

0157—  1 — -01 

0158 — 0 00 

0159 — 0 00 

0 1 68*  -XT 0 23 

0161—  + 33 


102 


104 


■iaag m 


1 


I 


I 


II. 

I 


0972— 

0973— - 

0974 — 

0975— - 

0976 — 

0977— 

0978— 

0979  — 

0980  — 

0981  — 

0982 — 

0983 — 
0984  — 
9985  — 
0986— 
0987  — 
C988— 

0989- 

0990- 

0991- 

0992- 

0993- 
9994- 

0995- 

0996- 

0997- 

0998- 
3999- 
1000- 
1001- 

1002- 

1003- 

1004- 

1005- 
1 006- 

1007- 

1008- 
1009- 
1010“ 
1011“ 
1012" 

1013- 

1014- 

1015- 

1016- 

1017 

1 0 1 8 

1019 

1020 
1021 
1022 
1 023 

1024 

1025 


UP — 

C.  1 

1026—  UP — 

XFR— 

6 7 

1 027 — 2 — 

-07 

1028— DIV— 

o 

&_ 

-02 

1029—  1 — 

*£<  Y — 

-52 

1038— XEY— 

CUT  — 

-47 

1031—  - — 

$/R— 

—1 

“ 4* 

1032 — DN— 

cut— 

-47 

1033— CHS— 

•CUT— 

-47 

1034—  X -• 

•XFR— 

-6  < 

1035— RUP— 

I till— 

-31 

1 036- -XEY-' 

b — 

-14 

1037— XFR- 

-XT0 — 

-23 

1038  1 — 

. 2 — 

-02 

1039—  2 - 

-CLX — 

-37 

1040 — XFR- 

- 1 -- 

-01 

1041—  - - 

-XTQ — 

-23 

1042--  2 - 

- + — 

— 3 

1043—  X - 

- b — 

-14 

1044 — I'H- 

-XFR— 

-67 

1045— XFR- 

-IND-- 

-31 

1046 — + - 

- b — 

-14 

1047—  1 " 

-XTQ — 

— 2 "-S 

C—  V? 

1048 — D IV- 

_ i -- 

-01 

1049— YTO- 

- UP — 

-27 

1050—  6 - 

-XTQ — 

-02 

-23 

1851—  9 - 
1052— XFR- 

- + - 

-33 

1053 — 3 - 

- b - 

-14 

1054— XFR- 

-XFR- 

- fu7 

1055 

-IHD- 

-31 

1056—  2 - 

- b - 

— 14 

1057—  r - 

-XT0- 

--23 

1058—  UP- 

- 0 - 

— 88 

1059—  5 - 

—34 

106Q — X - 

- 1 “ 

— Ol 

1061 — I'H- 

-XTO- 

—23 

1062 — XTQ- 

--34 

1063—  7 - 

- b - 

— 14 

1064 — 5 - 

-XFR- 

— 67 

1065— XFR- 

■j  — 

t- 

--02 

1066—  6 

- UP- 

— 27 

1067 — 9 

-XFR- 

— 67 

1068—  UP- 

- 3 - 

— 03 

1 069 — XFR 

^ . 

— 34 

1070—  7 

— DH- 

— 25 

1071—  5 

— DIV- 

— 35 

1072-'-  X 

— YTO- 

—40 

1073— YTO 

--07 

1074—  6 

--  8.  - 

—03 

1075—  1 

— XFR- 

--67 

1076 — XTQ 

-j 

— 07 

1077—  + 

4 - 

— 04 

1078—  7 

— XSQ- 

— 12 

1079—  1 

27 

02 

35 

01 

30 


-34 

-25 

-32 

-36 

-22 

-SO 

-67 

-07 

-02 

-67 

-34 

-02 

-36 

-25 

_ c 

O i 

-33 
-01 
-35 
-40 
-06 
-11 
-67 
-03 
—67 
-34 
—02 
—76 
—27 
—05 
— 36 


— 0 r 


■67 

•06 

■11 


— t* » 

—07 

— 05 
--36 
—40 


•06 

•01 


—33 
— 07 
— 0 1 


1 080  — X' 36 

1881—  UP 27 

1082—  XFR— 67 

1083—  7 07 

1084—  4 04 

1085—  X 36 

1 086—  RUP 22 

1087  — XEY--- 30 
1088—  2 02 

1089 —  DIV 35 

1090—  DN 25 

1091—  + —33 

1 092-- VTO 40 

1093 — + 33 

1094—  7 —07 

1095—  2 02 

1096 —  XFR 67 

109  7 — 6 06 

1098 — 1 01 

1099  — XTQ 23 

1100 — + 33 

HOI—  7 —07 

1102—  4 04 

1 103—  XFR— 67 

1104—  5 05 

1105—  9 —11 

1106—  UP 27 

1107—  CLX— 37 

1108—  X<Y— 52 

1109—  1 81 

1110—  1 01 

1111—  2 02 

1112—  4 —04 

1113—  XFR— 67 

1114—  7 —0? 
111  5 — 2 —02 

1 116—  CHS— -32 

1117—  UP— 27 

1118—  XFR--- 67 

1119—  7 —07 

1120—  1 01 

1121—  FMT  — 42 

1122—  1 —01 

1123—  DN— -25 
1 124  — XFR 67 

1125—  2 —02 

1126—  UP 27 

1 127— XFR— 67 

1128—  7 —07 

1129—  2 02 

1 ISO  — X<Y 52 

1131—  0 GO 

1 1 32--  9 —1  1 
1133—  3 --03 


L 

^ — ..I.  


r 


r 


1296—  M - 

—70 

1350- 

-EEX-- 

-26 

1 4 0 4 - ■ 

- 3 

-03 

1297—  ft  - 

-62 

1351- 

- 4 __ 

-04 

1405- 

- 0 -- 

-00 

1298"  YE- 

--24 

1352- 

-FMT  — 

-42 

1406- 

- 0 -- 

-00 

1299— CHT- 

-47 

1353- 

- 1 — 

-01 

14  07- 

- UP  — 

-27 

1300 — o.  - 

— 13 

1354- 

- tr 

•— 

-05 

1408- 

- ft.  — 

-13 

1301—  H - 

— 62 

1355- 

- 5 -- 

-05 

1409- 

-FHT-- 

-42 

1302 — XFR- 

—67 

1356- 

- 0 

-OO 

1410- 

- 1 -w 

-01 

1303 — CUT- 

--47 

1357- 

- 0 — 

-GO 

1411- 

- UP— 

-27 

1304—  ft  - 

— 62 

1358- 

-FMT  — 

-42 

1412- 

- b — 

-14 

1305—  N - 

— 73 

1359- 

- 1 — 

-01 

1413- 

-FMT— 

-42 

1306 — i j - 

— 15 

1360- 

- 6 — 

-06 

1414- 

- 1 — 

-01 

1 307--FMT- 

— 42 

1361- 

” a 

-21 

1415- 

-PUT  — 

-45 

1308--STP- 

__41 

1362- 

- 1 -- 

-01 

1416- 

-A  — 

-06 

1309 — PMT- 

--45 

1363- 

- 0 -- 

-00 

1417- 

-EEX-- 

-26 

1310— XFR- 

--67 

1364- 

-FMT  — 

-42 

1418- 

- 4 -- 

-04 

1311—  X - 

— 36 

1365- 

- 1 -- 

-01 

1419- 

-XTO  — 

“ C.O 

1312—  ft  - 

— 13 

1366- 

- 1 -- 

-01 

1420- 

- + — 

-33 

1313 — XTO- 

— 23 

1367- 

- 1 — 

-01 
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Transmission  Loss  / Probability  of  Detection  Program 
flow  charts  and  sample  outputs 
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1255— -  0 00 

1256 — 4 04 

1257—  XTO— 23 

1258 — 6 00 

1259 — 8 10 

1260 — 8 10 

1261—  1 —01 

1262—  . —21 

1263 — 2 02 

1264 — 8 10 

1265—  XTO  — 23 

1266 — 0 OO 

1267—  8 —10 

1268—  9 —11 

1 2 69--  1 01 

1 270- -  . —21 

1271—  6 —“06 

1 2 7 2 - - 4 04 

1273  — XTO 23 

1 274- -  9 —1  1 

1275 — 0 0O 

1276  — CUT— -47 

1277--  2 02 

1278  — EEX  — -26 

1279— -  3 03 

1280—  XTO 23 

1281—  9 “—11 

1 2 y 2 — 6 u 6 

1283— XTO— -23 
1284  — O OO 

1285—  9 11 

1286 — 7 07 

1287—  1 —01 

1 288—  XTO 23 

1289— -  3 03 

1 290—  XT  O 23 

1291— -  9 11 

1292 — 8 10 

1293—  XTO 23 

1294 — 8 IP 

1295 — O OO 
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1 732-- 

X -- 

-36 

1836- 

»*i 

IJ 

-10 

1 890- 

-XFR— 

- 

1 783-- 

YTO— 

-48 

1837- 

-XFR— 

-67 

1891- 

_ "j? 

- 

1784  — 

- -- 

-34 

1838- 

__  “7 

»‘ 

-OF 

1892- 

- 4 — 

- 

1 785-- 

pi 

-07 

1.839- 

-7 

c 

-07 

1893- 

- - — 

- 

1786  — 

5 — 

-05 

1840- 

-GTO — 

-44 

1894- 

-YTO— 

- 

1 787-- 

GHT — 

-47 

1841- 

-S/R  — 

“7  *"* 
“ t 1* 

i 895- 

— 7 ’ 

- 

1788— 

1 -- 

-01 

1842- 

- 1 — 

-01 

1896- 

_ 6 — 

- 

1783  — 

XTO— 

-23 

1843- 

- 0 

-00 

1897- 

-XFR— 

- 

1790  — 

+ --- 

-33 

1844- 

-XFR-- 

-67 

1 898- 

- 6 -- 

- 

1791  — 

9 — 

-11 

1845- 

- 7 -- 

-07 

1899- 

<L  ” 

- 

1 7 9 2 - - 

-10 

1846- 

-3  — 

-03 

19O0- 

-D IV  — 

- 

1 793-- 

GTO — 

-44 

1847- 

-GTO — 

-44 

1901- 

- G — 

- 

1 7 9 4 - - 

1 -- 

-01 

1848- 

-S/R— 

— 77 

i i 

1902- 

-YTO  — 

1795- 

r% 

»_•  ” ” 

-10 

1849- 

- 1 -- 

-01 

1903- 

- 0 - - 

- 

1796- 

2 — 

-02 

1850- 

- 0 -- 

- 0 O 

1904- 

-CLX-- 

- 

1797- 

£. 

-02 

1851- 

- mi— 

-25 

1905- 

-X*Y— 

- 

1798- 

1 — 

-01 

1852- 

-4-  — 

-33 

1986- 

- 1 -- 

- 

1799- 

X TO-- 

-23 

1853- 

-XFR— 

_67 

1907- 

— Q 

•* 

- 

1000- 

9 -- 

-11 

1854- 

- 7 -- 

-07 

1908- 

- 8 - - 

- 

1801- 

i"i 

ij 

-10 

1855- 

- 8 - - 

-10 

1909- 

_ r« 

- 

1802- 

XTO— 

-23 

1856- 

-XFR— 

-67 

1910- 

- q — 

m" 

- 

1883- 

O 

-03 

1857- 

— ,f  _ 
i’t 

-36 

1911- 

-5  — 

- 

1804- 

CMT-- 

-47 

1858- 

- 7 -- 

-07 

1912- 

— v(  J | J 

- 

1805- 

XFR  — 

- 6 7 

1859- 

__  -7 

i 

-07 

1913- 

- 0.  — 

- 

1886- 

9 — 

-11 

I860- 

-33 

1914- 

IJ  “ 

- 

1807- 

7 -- 

-07 

1861- 

-YTO— 

-40 

1915- 

-9  — 

- 

1808- 

X TO-- 

-23 

1862- 

- 7 — 

-07 

1916- 

-XTO  — 

- 

1809- 

9 -- 

-11 

1863- 

- 5 — 

-05 

1917- 

- k -- 

•- 

1810- 

6 — 

-06 

1864- 

-GTO  — 

-44 

1918- 

-XFR — 

- 

1811- 

XFR — 

-67 

1865- 

-S/R— 

_ “7  "7 
t 1 

1919- 

- 0 - - 

- 

1812- 

- S --- 

-10 

1866- 

- 1 -- 

-01 

1920- 

- UP-- 

- 

1813- 

-XTO— 

_ 2 O 

1867- 

- 7 — 

-07 

1921- 

-XFR— 

- 

1814- 

- + — 

-33 

1 8 6 8 — 

-XFR— 

-67 

1922- 

- 9 -- 

- 

1815- 

- Q 

•r 

-11 

1869- 

-07 

1923- 

- 0 - - 

- 

1816- 

«.  o 

-03 

1870- 

- 4 -- 

-04 

1924- 

- X > Y - - 

~ 

1817- 

-XFR — 

-67 

1871- 

- U H S — 

-32 

1925- 

- 1 -- 

- 

1818- 

-9  — 

-11 

1872- 

- UP— 

-27 

1926- 

- 9 -- 

- 

1819- 

-XTO  — 

— 0*7* 

Cm  V 

1873- 

-XFR— 

-67 

1927- 

- 3 — 

-- 

1828- 

- -r  -- 

_ 3 ?. 

1874- 

- 7 -- 

-07 

1928- 

i*. 

- £•  ““ 

- 

1 82 1 - 

- 0 -- 

-ii 

i r»  cr 
1 O ( J - 

- 7 -- 

-07 

1929- 

-EEX-- 

- 

1822- 

- 4 - - 

-04 

1876- 

-FMT  — 

-42 

1930- 

- 2 — 

- 

1823- 

-GTO— 

-44 

1877- 

- 1 -- 

-01 

1931- 

••XTO— 

- 

1824- 

-S/R-- 

-77 

1878- 

- IJP— 

-27 

1932- 

- a - - 

- 

1825- 

- 1 -- 

-01 

1879- 

-FMT— 

-42 

1933- 

-GTO  — 

- 

1826- 

•1 

- f — 

-07 

1888- 

--  1 -- 

-01 

1934- 

- 1 -- 

- 

1827- 

-CL  2 — 

_ V 

I 

1881- 

- 4 -- 

-04 

1935- 

- 9 — 

- 

1828- 

- IJp— 

-27 

1882- 

-5  — 

-05 

1936-- 

- 7 -- 

- 

1829- 

-XFR — 

- 6 7 

1883- 

- 0 - - 

-00 

1937- 

- 1 -- 

- 

1830- 

-02 

1884- 

-XT0-- 

~ o •:* 

Cm  *-* 

1938- 

-XFR  — 

- 

1831- 

- tr  _ 

-05 

1 OC«tr 

- o.  - - 

-13 

1939- 

-I  HD— 

- 

1832- 

• ,»  _ » » 

1 

-50 

1 o »:•  u _ 
X 'j  C' 

-XFR— 

-67 

1940- 

- b — 

- 

1 833- 

- l — 

-01 

1807- 

-06 

1941  - 

- UP-- 

- 

1834- 

„ i*«  m . , 

-10 

1888- 

_ i 

A 

-01 

1 942- 

-XFR— 

- 

1835- 

- I > — -* 

-06 

1889- 

~ UP  — 

~ j 7 

u 1 

1 943- 

- 0 - - 

- 

t*  f 

07 

04 

34 
40 
07 
06 
67 
06 
02 

35 
15 
40 
00 

i- 

50 
01 
11 
10 
10 
1 1 

05 

■j  o 

13 

10 

1 1 

00 

n 

K *7 

t*  I 

00 

■‘j7 

i_  i 

67 

1 1 

00 

cr 

._i , 

01 
11 
03 
10 
26 
02 


‘t*T 

01 

1 1 

07 

01 

"i 

31 
14 
27 
6 7 
00 


PEAPS 

FORTRAN  Program  Listing 
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c 

c 

c 


PEAPS  (FORTRAN  VERSION) 


DIMENSION 

DIMENSION 

DIMENSION 


HEAD( 20) , 
SGRAD( 15) , 
I H(  6 ) 


D(  15) , 
CD(  6 ) 


INTEGER  * 2 IBUF(6l,  121) 

DATA  X/. 12, .25, .38, . 52, .67, 
INTEGER  * 2 IB/*  */ 

INTEGER  * 2 IDOT/'.  '/ 

INTEGER  * 2 1ST AR/ * * •/ 


T(  15) , SALIN(  15)  , S( 15) 
, CS ( 6 ) , PR(  13) , X( 8 ) 

.84, 1.06,1.28/ 


C 

C 

C 


C 

C 

C 


C 

c 

c 


PI  = 3.  141592 
RAD  = PI  / 180.0 

90  DO  91  1=1,6 
CD ( I } = 0.0 

91  csm  = o.o 

READ  AND  PRINT  HEADING  CARD  (80  COLUMNS  MAX) 

READ  (5,  900,  END  = 9999)  HEAD 

900  FORMAT  (20A4) 

WRITE  (6,  901)  HEAD 

901  FORMAT  (1H1,  20A4) 

READ  AND  PRIMT  INPUT  DATA  CARDS 

READ  (5,  902)  NPTS,  MODE,  IBT , IWH,  ISS,  IFOM,  SD,  RD, 

❖ SF 

902  FORMAT  (615,  3F10.2) 

FOM  = FLOAT (IFOM) 

SS  = FLOAT (ISS) 

WRITE  (6,  903)  NPTS,  MODE,  IBT,  IWH,  ISS,  IFOM,  SD, 

* RD,  SF 
( 1H0,  22HNUM3ER  OF  SSP  POINTS  =,  15// 

19H  RAY  PLOTTER  MODE  =,  15  // 14.1  BOTTOM  TYPE  =, 

15  //14H  WAVE  HEIGHT  =,  15,  3H  FI  // 

15H  SYSTEM  SIGMA  =,  15  //6H  FOM  =,  1 5 // 

13H  SOURCE  DEPTH  , F10.2,  3H  FT  // 

15H  RECEIVER  DEPTH  , F10.2,  3H  FT  // 

17H  SOURCE  FREQUENCY  , F10.2,  3H  HZ  ) 

(6,  901)  HEAD 
(6,  904) 

( 1 HO » 5X , 5HP0I NT , IX , 5HDE PTH ,4X , 1 lHTEMPERATURt • 
IX,  8HSALINITY,  IX,  11KS0UND  SPEED,  IX, 

8HGRADI ENT ) 


903  FORMAT 
* 

* 

* 

* 

WRITE 

WRITE 

904  FORMAT 
* 

* 


SPEED  OF  SOUND  PROGRAM 


905 

996 


Id), 


$ ALI N ( 


DO  102  I = 1,  NPTS 
READ  (5,  905)  0(1), 

FORMAT  (3F10.2) 

WRITE  ( 6,  906)  I , D ( I ) , T ( I ) , 
FORMAT  ( 1H0,  110,  3F10.2) 

T ( I ) = ( T ( I ) - 32.0)  - 5.0  / 9.0 
= 0(1)  * 0.3048 
0(1)  * D ( I ) 

D2  * D2 
T(  I)  * T ( I ) 

T2  * T( I) 

SALIN(I)  - 35.0 
SI  * SI 

1449.3  + 1.5848E-2  * D(I) 


SALlNt I) 


D ( 1 ) 
D2  = 
D4  = 
T 2 = 
T3  = 

51  = 

52  = 
S(I) 


* 

* 

* 


907 


S ( I ) = 
D ( 1 ) * 
WRITE 


- 3.46E-16  * D4 
+ 2.604E-4  * T3 
+ 1.534E-4  * T 2 

- 1.2E-2  * SI  * 
S(I)  / 0.3048 
D( I ) / 0.3048 

(6,  907)  S(  I ) 


+ 4.587  * 
<■  1.19* 
* D ( I ) - 
T ( 2 ) 


+ 1.572E-7  * D2 
Til)  - 5.356E-2  * T2 
SI  + 9.6E-2  * S2 
7.19E-9  * T ( I # * D 2 


FORMAT  (1H-,  40 X , F10.2) 
IF  (I  .GT.  1)  30  TO  100 
SGRADI 1 ) =0.0 
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ooo  ooo 


GO  TO  101 

100  J = I - 1 

SGRAD( I)  = ( S< I)  - S(J)J 

101  WRITE  (6,  908)  SGRAO(I) 
908  FORMAT  (1H0,  50X,  F10.4) 

102  CONTINUE 


/ (0(1)  - DC  J ) J 


SPEED  GF  SOUND  PROFILE  CRITICAL  POINTS  PROGRAM 


I = 0 
I FLAG  = 0 
I RUN  = -1 
J = 1 

K = 2 

NG  = 0 
NTEST  = -1 

103  I = I 4-  1 

J = J + 1 

NTEST  = NTEST 

104  I RUM  = I RUN  4- 
IF  (J  .GT.  NPTS ) GO  TO 
IF  (SI  I ) .EQ.  S( J)  ) GO 
IF  (NTEST  .EQ.  0)  GO  TO 


4- 

1 


NTEST  = 

cs(k)  = sin 

CD(K)  = 0(11 

K = K 4- 

NG  = NG 


113 

TO  103 
105 


1 

4-  1 


IF  ( S ( J ) .LT.  S( I ))  GO 
I FLAG  = 1 

105  IF  (IFLAG  .EQ.  1)  GO  TO 

106  IF  (S(J  ) .LT.  S( I)  ) GO 
IF  ( IRUN  .EQ.  0) 

CS( K)  = S 1 1 ) 

CD(K)  = 0(1) 

1 


TO  105 


GO  TO 


110 

TO  108 
107 


4- 

1 


GO  TO  109 


4-  1 


K = K 4- 

NG  = NG 
IFLAG  = 

GO  TO  104 

107  IFLAG  = 1 
GO  TO  104 

108  IF  (IRUN  .GT.  0) 

CS( K)  = S(I) 

CD  ( K)  =0(1) 

K = K 4-  1 

NG  = NG  4-  1 

109  IRUN  = IRUN 

1 = 14-1 
J = J 4-  1 

IF  (J  . GT.  NPTS)  GO  TO 
GO  TO  106 

110  IFLAG  = 0 

IF  (S(J)  .GT.  S ( I ) ) GO 
CS(K)  = S ( I ) 

C0( K)  = D ( I ) 

K = K 4-  1 

NG  = NG  4- 
GO  TO  104 

111  I = I 4-  1 
J = J 4-  1 

IRUN  = IRUN  f 1 

IF  (J  .LE.  NPTS)  GO  TO 

112  CS(K)  = S( I; 

CO(K)  = D ( I) 

4-  1 


112 

TO  111 


1 


110 


NG  = NG 

113  CS<  1)  = S(  1) 

CD( 1)  = D( 1) 

CS( 6)  = S( NPTS ) 
CD ( 6 ) = D(MPTS) 


SCUND  SPEED  PROFILE  PROGRAM 
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WRITE  <6,  901)  HEAD 

909  FORMAT1 (IhO?935HSOUND  SPEED  PROFILE  CRITICAL  POINTS  // 

* 5X,  6H  DEPTH » AX,  14HSPEE0  OF  SOUND  ) 

DO  114  I = It  6 

114  WRITE  (6,  910)  C0(I),  CS(I) 

910  FORMAT  (1H0,  2F10.2) 

DO  115  I = 1,  61 

DO  115  J = It  121 

115  I BUF( I , J)  = 18 
WRITE  (6,  901)  HEAD 

911  FORMAT*  'll  HO  ^37HSPE  ED  OF  SOUND  (FT/SEC)  VS  DEPTH  (FT)) 
SMIN  = CS(1) 

SMAX  = CS ( 6) 

DO  116  I = 1,  6 

IF  (CS(I)  .EG.  0.0)  GO  TO  116 

IF  (CS(I)  .LT.  SMIN)  SMIN  = CS ( I ) 

IF  tCS(  1 1 .GT.  SMAX)  SMAX  = CS ( I ) 

116  CONTINUE  „ „ 

PRI  = ( SMAX  - SMIN)  / 12.0 
PR( 1 ) = SMIN 

V = PR ( 1 ) 

DO  117  I = 2,  13 

V = V + PRI 

117  PR(I1  = V 

WRITE  (6,  912)  ( PR( I ) , I = 1,  13) 

912  FORMAT  (F8.1,  12F10.1) 

DG  121  I = 1,  61  _ ^ 

IF  (MOO  ((1-1),  10)  .EQ.  0}  GO  TO  119 
DO  118  J = It  121,  10 

118  I BUF(  I,  J)  = IDOT 
GO  TO  121 

119  DO  120  J = 1,  121 

120  I BUF( I,  J)  = IDOT 

121  CONTINUE 
WRITE  (6,  913)  ( I B:JF  ( 3 , 

913  FORMAT  ( 1H  , 5X,  121A1) 

PRI  = ( PR ( lb)  - PR  till  / 

/ 60.0 


J) , J = 1,  121) 
1 2 U r.  0 


PPDI  = D(NPTS) 

J = NPTS  - 1 
DO  123  I = 1,  J 
PPD  = 0(1) 

PPR  = S ( I ) 

PPR I = SGRADU  * 1)  * PPDI 

122  I PPD  = (PPD  / PPDI)  * 1*5 
I PPR  = (PPR  - PR( 1 ) ) / PRI 
IBUFdPPD,  IPPR ) = ISTAP. 

PPD  = PPD  + PPDI 

IF  (PPD  .GT.  D(  I <-  15)  GO  TO  123 
PPR  - PPR  + PPR I 
GO  TO  122 

123  CONTINUE 
MND  = O(NPTS)  / 6.0 
MNDP  = 0 
DO  126  I = It  61 
IF  ( I . EQ.  1 ) GO  TO  125 
IF  (MOD ( ( I~1 ) , 10)  .EQ.  0 GO  i0,124 
WRITE  (6,  913)  ( I BUF(  I , J ) , J ~- 
GO  TO  126 

124  MNDP  = MNDP  *■  MND 

125  WRITE  (6,  914)  MNDP,  (IBUF(1,J), 

914  FORMAT  ( IH  , 15,  121A1 ) 

126  CONTINUE 

RAY  TRACE  PROGRAM 
I FLAG  =0 

IF  (MODE  . NE . 2)  GO  TO  121 

READ  IN  MOOE  2 RAY  TRACE  EMISSION  ANGLES 


121) 


J = 1,  121) 
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READ  ( 5 » 905)  THM IN*  DTH,  THMAX 

127  IF  (MODE  ,EQ.  1)  GO  TO  134 
DO  128  I = It  61 

DO  128  J = 1,  121 

128  IBUF  ( I , J)  = 18 
WRITE  (6;  901)  HEAD 
WRITE  (6,  915) 

915  FORMAT  (iHO,  24HRANGE  (NM)  VS  DEPTH  (FT)  ) 

PRI  = 5.0 

V = -5.0 

DO  129  I = 1,  13 

V = V + PRI 

129  PRC  I ) = V 

WRITE  (6,  912)  (PR(I)t  I = It  133 
DO  133  I = lt  61 

IF  (MOD  ((I  - ‘ 1 ) , 10)  .EQ.  0)  GO  TO  131 
DO  130  J = 1,  121,  10 

130  IBUF( I,  J)  = I DOT 
GO  TO  133 

131  DO  132  J = 1,  121 

132  IBUF( I,  J)  = IDOT 

133  CONTINUE 

PRI  = 3000.0 
RMAX  = 360000.0 

WRITE  (6,  913)  C I8UFC  3 , J),  J = 1,  121) 

CONVERGENCE  ZONE  TEST 

IF  (MODE  .EQ.  2)  GO  TO  219 

134  IF  (CS ( 6 ) . GT.  C S ( 2} 1 GO  TO  136 

135  WRITE  (6,  916) 

916  FORMAT  (IHO,  19HN0  CONVERGENCE  ZONE  ) 

ICZ  = 0 

GO  TO  139 

136  M = 2 

137  M = M + i 

IF  (CS(M)  .EQ.  0.0)  GO  TO  135 
IF  (CS( 2 ) .GT,  C S ( M ) ) GO  TO  138 
IF  (M  .EQ.  6)  GO  TO  135 
GO  TO  137 

138  ICZ  = 1 

139  M = 0 

200  M = M + 1 

IF  ( SD  .GE.  D ( M ) ) GO  TO  200 

SV  = ((SD  - D(M  - 1))  * SGRAD(M))  + S(M  - 1) 

THP  = (2.0  - (2.0  - SV  / CS ( 2 ) ) ) **  0.5 
THPP  = (A8S  (2.0  - (2.0  * SV  / CS(6)))>  **  0.5 
IF  (ICZ  .EQ.  0)  GO  TO  227 
THMIN  = THPP 
I PLOT  = 1 
I SR  = 1 

201  IF  (MODE  . EQ.  1)  I PLOT  = 0 
TH  = THMIN 

Z = SD 
R = 0.0 

202  M = 0 

203  M = M + 1 

IF  (DIM)  .GT,  Z)  GO  TO  204 
IF  (M  .EQ.  NPTS)  GO  TO  208 
GO  TO  203 

204  IF  (Z  .LT,  0.0)  GO  TO  208 

205  RADIUS  = (-SGRAD(M)  * (1.0  - ( ( TH  * TH)  / 2.0)))  / 

* ((SGRAD(M)  * (Z  - D(M  - 1) ) + S(M  - 1))) 

DR  = 5.0  =?  ( ( D ( M ) - D ( M - l>)  **  0.5) 

206  RDR  = DR  * RADIUS 
R = R * DR 

Z = Z + ((RADIUS  * (DR  * DR))  / 2.0)  + ( TH  * DR) 

TH  = TH  + RDR 

IF  (I  PLOT  .EQ.  0)  GO  TO  207 
I PPD  = ( /.  / PPOl  ; <•  1.5 
IF  ( IPPD  .oT.  61 ) GO  TO  207 
I PPR  = R / PRI 
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. GT . 


207 


208 

209 


210 


211 


212 


213 


214 

215 


216 


217 


218 


219 

220 


221 

222 


223 


IF  { I PPR 
I3UFUPPD,  I PPR ) 
IF  (Z  .LT.  DIM  - 
< D<  M J .LT.  Z) 

( RMAX  .GE.  R) 
TO  {209,  212, 
(C(NPTS)  .LT. 

( IFLAG  .EQ.  i) 


IF 

T 


202 


Ir 

r, : 


GO 


IF 

IF 


THMIN 
I PLOT 
I SR  = 
GO  TO 
THMIN 
1 PLOT 
I FLAG 
I SR  = 
GO  TO 
THBB  = 


121)  GO  TO  207 
= ISTAR 
D)  GO  TO 
GO  TO  202 
GO  TO  206 
213,  216,  217, 
Z)  GO  TO  210 
GO  TO  211 


222,  220) 


IFLAG 
ICZ  = 
! PLOT 
CZMAX 
I SR  = 
GO  TO 
THMIN 
I PLOT 
ISR  = 
GO  TO 
IF  (R 
THMIN 
I PLOT 
ISR  - 
GO  TO 
IF  (R 
CZM1N 
THMIN 
’PLOT 
ISR  = 
GG  TO 
IF  (R 
THMIN 
I PLOT 
ISR  = 
GO  TO 
IF  (R 
CZMIN 
THMIN 
ISR  = 
GO  TO 
THMIN 
I PLOT 
ISR 


= THMIN  + (0.5  * RAD) 
= 1 
1 

201 

= THMIN  - (0.1  * RAD) 
= 0 
= 1 
1 

201 

THMIN 


CZMAX)  CZMAX  = R 


= 0 
1 

= 1 
= R 
2 

201 
= TnP 
= 1 
3 

201 

.GT.  (CZMAX  / 2.0))  GO  TO  214 
= THMIN  + (0.1  * RAD) 

= 0 

3 

201 
.GT, 

= R 

= THMIN 
= 1 

4 

201 

.LT.  CZMIN)  GO 
= THMIN  - (0.4 
= 0 

5 

201 

.GT.  CZMIN)  GO  TO  218 
= R 

= THMIN  <-  (0.1 
5 

201 

= THMIN  - (0.1 
= 1 
t> 

201 

.EQ.  1 ) GO  TO  222 


+ (0.5 


RAD) 


TO  215 
* RAD) 


* RAD) 


* RAD) 


GO  TO  ___ 

IF  (IFLAG 

I C|,'  - f 

GO  TO  201 
THMIN  = THMIN 
IF  (THMIN  .LC 
ISR  = 1 
GO  TO  223 
MODE  = 1 
GO  TO  L 34 
1 

= CZMIN 
= CZMAX 
= l 
2 

= 0 
I = i 


+ DTH 
, THMAX ) 


GO  TO  219 


IBB  = 
CZMIN 
CZMAX 
I PLOT 
ISR  = 
MNDP 
DO  226 


/ 

/ 


3.0 

3.0 


IF  (I  . EQ.  ll 
.IF  (MUDCU-i), 
WRITE  (6,,  913) 


61 

GO  TO  225 
10)  .HU. 
( 1 3UF(  I , 


0) 
J)  , 


GO  TO  224 
J " 1,  121) 


ISR 


.158 


GO  TO  226 

224  MNOP  = MNOP  4 MND 

225  WRITE  (6,  914)  MNOP, 

226  CONTINUE 

GO  TO  ( 221,  227) , ISR 

227  CONTINUE 


( I6UF( I,J) , 


121  ) 


TRANSMISSION  LOSS  PROGRAM 
B = FOM 

IF(80.GE. IFOM ) B = 80.0 
IF{ 100. LE • IFOM)  B = 100.0 
B = (B  4 15.0)  / 20.0 
B = 10.0**B 
IR  = B/20000.0 
RMAX  = (IR  4 1)  * 20000 
WRITE  16,  901)  HEAD 
WRITE  (6,  917) 

FORMAT  (1H0,  10HRANGE  (NM) 
DO  228  I = 1,  61 
DO  228  J = 1,  121 
IBUF(I,  J)  = IB 
PRI  = RMAX  / 24000.0 


4HPDET) 


= I DOT 
IPFCM) 


1,  121 
= I DOT 
IPFOM) 


.EQ.  0)  GO  TO  231 
10 

IBUFd,  J)  = ISTAR 


IBUF(I,  J)  = ISTAR 


120.0 


WRITE  16,  901)  HEAD 
WRITE  (6,  917) 

917  FORMAT  (1H0,  10HRANGE  (NM)  , IX,  5HTL0SS  , 4X,  4HPDET) 
DO  228  I = 1,  61 
DO  228  J = 1,  121 

228  IBUFd,  J)  = IB 

PRI  = RMAX  / 24000.0 
PR ( 1 ) = 0.0 

V = 0.0 

DO  229  I = 2,  13 

V = V 4 PRI 

229  PR ( I ) = V 

IPFOM  = ( ( IFOM  - 40)  / 2) 

DO  233  I = 1,  61 

IF  (MOD(d-l),  10)  .EQ.  0)  GO  TO  231 
DO  230  J = 1,  121,  10 
I BUF( I,  J)  = I DOT 

IF  (I  .EQ.  IPFCM)  IBUFd,  J)  = ISTAR 

230  CONTINUE 
GO  TO  233 

231  DO  232  J = 1,  121 
IBUFd,  J)  = I DOT 

IF  (I  .EQ.  IPFOM)  IBUF(I,  J)  = ISTAR 

232  CONTINUE 

233  CONTINUE 

PRI  = RMAX  / 120.0 
PPDI  = 2.0 

I F (CS ( 2 ) .IT.CSt 1 ) ) GO  TO  304 
IF ( SO. GT.CDt  2 ) ) GO  TO  303 
I F( KD.GT.CDl 2) ) GO  TO  301 
ISD  = 1 
IRAP  = 0 

RS  = ARC0S(CS(l)/CS(2))*(2.0/3.0)*(CS(l)*(CD(2)-CDd)) 

* ) / (CS  ( 2)  - CS  ( 1)  ) 

RT  = CD ( 2 ) /( 6 .O^ARCOS ( SV/CS ( 2) ) ) 

SSL  = (0.8/RS)*( (WH*SF/1000.0)**0.71) 

G = AB  S ( { CS ( 31 -C  S ( 2 > ) / ( CD ( 3 ) -CD  £ 2 ) ) ) 

DFL  = 744. 0*( (SF/1000.0)**(~5. 0/3.0) )*(G**(-1. 0/3.0) )* 

* ( CD( 2)  **  (-3.0)) 

GO  TO  306 

301  IRAP  = 1 

RL2  = CS(2)*(CD(3)-CD(2) ) / ( CS( 2) -CS( 3) ) 

I F{ SD.GT ,CD( 2 ) ) GO  TO  302 
RL2  = RL2*(RD-CD(2))/4.5 
GO  TO  305 

302  RL2  = RL2*(SD-CD(2) )/4.5 
GO  TO  305 

303  IF(RD.GT.CD(2) ) GO  TO  304 
GO  TO  301 

304  IRAP  = 0 

305  ISD  = 0 
GO  TO  307 

306  IF(NG.F.Q.l)  DFL  = 0.0 

307  1 F ( IC2.EQ.il  GO  TO  308 
GO  TO  310 

308  C2STRT  = CZMIN 
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C2ST0P  = CZMAX 
IFIIRAP.EQ.U  GO  TO  310 
I F t ISD.EQ.l ) GO  TO  310 
IDSC  = 1 
M = 1 

309  M = M + 1 

IF(CS(M+1).LT«CS(M) ) GO  TO  309 

RPT  = { CZMIN/  8.0)*( (CD(M)-CD(2) ) / < SO-CD  C 2 ) ))**0.5 
GO  TO  313 

310  FTEST  = 32.0  * (<SV  / SF)  **  0.7) 


I FI FTEST .GT. 

CD(  2 ) ) 

GO  TO  311 

GO  TO  312 

311 

ISO  = 0 

312 

I CSC  = 0 

313 

IFI6.LT. IBT) 

GO 

TO 

314 

IFI5.LT. I3T) 

GO 

TO 

315 

IFI4.LT. IBT) 

GO 

TO 

316 

IFI3.LT. IBT) 

GO 

TO 

317 

IFI2.LT.IBT) 
BT  = 0.66 
GO  TO  319 

GO 

TO 

318 

314  BT  = 2.0 
GO  TO  319 

315  BT  = 1.41 
GO  TO  319 

316  BT  = 1.18 
GO  TO  319 

317  BT  = 1.0 
GO  TO  319 

318  BT  = 0.85 

319  BLMA  = 30.0 

I F ( 3500 .O.LT.SF)  GO  TO  320 
I F ( 2000 .0. LT .SF ) GO  TO  321 
BLMA  = 50.0 

IFdOOO. O.LT.SF)  GO  TO  322 
IF(500. O.LT.SF)  GO  TO  323 
IF(LOO. O.LT.SF)  GO  TO  324 
B 1 = 0.2 
BO  = -0.83 
60  TO  325 

320  B1  **  0.35 
BO  = 7.12 
GO  TO  325 

321  B 1 = 0.28 
BO  = 8.21 
GO  TO  3 25 

322  B 1 = 0.19 
BO  = 2.09 
GO  TO  325 

323  B 1 = 0. 17 
60  = 1.97 
GO  TO  325 

324  B 1 = 0.  18 
80  = 2.06 

325  SFKHZ  = SF/1000.0 
SFKHZ2  = SFKHZ*SFKHZ 

A = ( (0.1*SFKHZ2)/( 1.0+SFKHZ2) ) + < ( 40* SFKHZ2 ) / { 410 0.0 
* + S FKHZ2 ) ) 

A = A / 1000.0 
ICZRC  = 1 
H = 1.0 

DRWRK  = 2000.0 

ORPERM  = 2000.0 

IF  IICZ  .EQ.  1)  GO  TO  3251 

IF  { CS ( 6 ) .GT.  SV)  GO  TO  4250 

THP  = ARCOSISV  / CS ( 6 ) ) 

3250  THBB  = THP 

3251 


R = 0.0 

THBL  "PI/  2.0 
326  I F ( THBB .GT . THBL } 

R = R + DRWRK 
THBL  = ATAN(  (R*dC0(6  )-SO)/3.0)  )/ 


H = H + 1.0 
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* (H*( ( 2.0*CD(6 )-SD)/3.0) ) ) 

THBLA  = THBL  * (180.0  / PI) 

IF  (BLMA  .LT.  THBLA)  THBLA  = BLMA 
BL  = H * BT  * ((THBLA  * 61)  «•  BO) 

STL  = (20.0*AL0G10(R/C0S(THBL) } f- 

* ( A* ( R/COS ( THBL ) ) ) 4-  BL  - 6.0) 

IF( ISO.EQ.O)  GO  TO  330 

I SR  = 1 

TL  = (A  4-  DFL)  * R 
IF  (RT  .LT.  R)  GO  TO  327 
TL  = TL  4-  20.0*ALOG10(R) 

GO  TO  ^28 

32?  TL  = TL  + 10.0*AL0G10(R)  + 10.0*AL0G10(  RT ) 

328  IF  (RS  .GT.  R)  GO  TO  329 
TL  = TL  + (SSL  * R) 

329  F = (A8S( STL  - TL i J / (-10.0) 

AD  = 10.0  * ( ALOGIO ( 1.0  4-  (10.0  **  F))} 

IF1STL.LT. TL)  TL  = STL 
STL  = TL  - AD 

GO  TO  (330,  334,  338,  339),  ISR 

330  IF(IRAP.EQ.O)  GO  TO  334 
ISR  = 2 

RAPTES  = ( RL2  **  0.5)  4-  ( RS  / 2.0) 

IF(R.LT.RAPTES)  GO  TO  331 
I RAP  = 0 
GO  TO  330 

331  IF  (SD.LT.RD)  GO  TO  332 

TL  = 20.0*ALGG10(  (SD*SD  4-  R*R)**0.5)  4-  A*((SD*SD  + 

* R * R)  **  0.5) 

GO  TO  333 

332  TL  = 20.0* ALOGIO ( ( RD*RD  4-  R*R)**0.5)  4-  A*((RD*RD  + 

* R * R)  **  0.5) 

333  GO  TO  329 

334  IF(ICZ.EQ.O)  GO  TO  338 
ISR  = 3 

TL  = 20 . 0*AL0G1 0 ( R ) f A*R 
IFIR.LT .CZSTRT ) GO  TO  337 
IF(R.GT. CZSTOP)  GO  TO  336 
IF( ICZRC.EQ.2)  GO  TO  335 
DRWRK  = ( CZSTCP  - CZSTRT) /10.0 
R = CZSTRT  - DRWRK 
ICZRC  = 2 
GO  TO  326 

335  RANGE  = R - CZSTRT 

GCZMAX  = 15.5  - 0. 0005* (CZ STOP  - CZSTRT) 

GB  = 0.36  / (GCZMAX*GCZMAX) 

RP  = RANGED GB**. 5)  / (0. 1 8*(CZST0P  - CZSTRT)) 

GCZ  = (RP/GB)  * EX P(  — ( RP*RP ) / (2.0*GB)) 

TL  = TL  - GCZ 
GO  TO  337 

336  CZSTRT  = CZSTRT  4-  CZMIN 
CZSTOP  = CZSTOP  4-  CZMAX 
DRWRK  = DR PERM 

-ICZRC  = 1 

337  GO  TO  329 

338  !F( IDSC.EQ.O)  GO  TO  339 
ISR  = 4 

TL  = 10.0*AL0G10(R)  + 10.0*AL0G10(RPT)  4-  A*R 
GO  TO  329 

339  PPD  = STL 
PPR  = R 

I PPD  a ((PPD  - 40.0)  / PPDI  ) 4-  1.5 
IF  ( I PPD  ,GT.  61)  I PPD  = 61 
I PPR  = (PPR  / PRI  ) 4-  1.0 
IF  { I PPR  .GT.  121)  I PPR  = 121 
IF  (IPPO  .GT.  60)  I PPD  = 60 
IBUFUPPD,  I PPR)  = I STAR 
C 

C PROBABILITY  OF  DETECTION  PROGRAM 
C 

PD  = 0.5 

SE  = FOM  - STL 
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XX  = A8SC SE/SS ' 

IFCXX.EQ.O.O)  GO  TO  343 
PO  = 0.95 
I = 8 

IFd.64.LE.XXJ  GO  TO  341 

340  IFCXX.GT.XC  I H GO  TO  342 
PD  = PO  - 0.05 

1 = 1-1 

IF(I.EQ.O)  GO  TO  342 
GO  TO  340 

341  PD  = 1.0 

342  IF{ SE.GT.O.OJ  GO  TO  343 
PD  = 1.0  - PO 

343  PPD  = 61.0  - (PO  * 10.0) 

I PPD  = PPD  + 0.5 

IF  d PPD  .GT.  61)  I PPD  = 61 
IBUFUPPD,  IPPR ) = ISTAR 
RF  = R / 2000.0 
WRITE  (6,  913)  RP,  STL,  PD 

918  FORMAT  (1H0,  3F10.2) 

IF(KMAX.GT.R)  GO  TO  326 
WRITE  (6,  901)  HEAD 
WRITE  (6,  919) 

919  FORMAT  (1H0,  35HTL0SS  1 DB)  ABOVE  PDET  (PERCENT)  VS  , 

* 10HRANGE  ( NM)  ) 

WRITE  (6,  912)  { PRC  I ) t I = it  13) 

WRITE  (6,  913)  (I3UFC3,  J)  , J = 1,  121) 

MNDP  = 40 

MND  =20 

DO  347  I = 1,  61 

IF  (I  .EQ.  1 ) GO  TO  346 

IF  (MOD ( ( I — 1 ) * 10)  .EG.  0)  GO  TO  345 

IF  Cl  .GT.  50)  GO  TO  344 

WRITE  (6,  913)  (IBUt-U,  J),  J = 1,  121) 

GO  TO  347 

344  IF  (I  .GT.  52)  GO  TC  345 
MNOP  = 100 

MND  = -10 

345  MNOP  = MNOP  + MND 

346  WRITE  (6,  914)  MNOP,  (IBUFCI,J;,  J = 1,  121) 

347  CGNTINUE 
GO  TO  90 

9999  STOP 
END 


1C2 


.-i*. 


APPENDIX  E 


PEAPS  STORAGE  LOCATIONS  IN  THE  HP9810A  CALCULATOR 

A.  The  operator  enters  all  angles  in  degrees.  The  zero 
degrees  direction  for  all  angles  is  horizontal  to  the  right. 
Clockwise  angles  (down)  from  the  horizontal  are  to  be  entered 
as  positive  (+)  angles.  Counter-clockv?ise  angles  (up)  from 
the  horizontal  are  to  be  entered  as  negative  (-)  angles. 

B.  All  angles  are  stored  and  used  in  radians.  Entered  angles 
are  automatically  converted  from  degrees  to  radians  for 
storage  by  the  program: 

degrees  = jjjg-  * radians 

radians  = * degrees 

C.  Enter  all  temperatures  in  degrees  fahrenheit  (°F) . 
Temperatures  are  automatically  converted  to  degrees 
centigrade  (°C)  for  use  in  the  SSP  program: 

°C  = | (°F  - 32) 

°F  = 32  + (|  * °C) 


163 


D.  Enter  depths  in  feet  (ft) . Depths  are  stored  in  feet 
and  are  converted  to  the  required  units  as  needed  in  the 
programs : 

meters  ='  0.3048  * feet 
feet  = meters/0.3048 
yards  = feet/3.0 

nautical  miles  (NM)  = yards/2000  = meters/1852. 

E.  Ranges  are  in  feet  in  the  ray  trace  program  and  are 
automatically  converted  to  units  of  yards  for  use  in  the 
transmission  loss/probability  of  detection  program. 

F.  The  following  pages  of  program  forms  indicate  the 
parameter  storage  locations  used  in  the  PBAPS  model.  The 
terms  listed  are  defined  in  the  Table  of  Symbols  and 
throughout  the  body  of  the  text. 
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